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FOREWORD 


This  report  was  prepared  by  Mr.  Elijah  W.  Turner  and 
Mr.  John  T.  Riechers,  Aerospace  Engineers  in  the  Loads  and  Criteria 
Group,  Structural  Integrity  Branch,  Structures  Division  at 
Wright-Patterson  AFB,  Ohio.  This  effort  was  performed  to  generate 
pressure  and  temperature  data  for  life  cycle  testing  of  the  B-1B 
transparency. 

This  effort  was  conducted  in  support  of  Project  24010505,  B-1B 
TRANSPARENCY  TEST,  which  is  managed  by  Vehicle  Subsystems  Division 
(WL/FIV) ,  Mr.  Richard  Smith,  Project  Engineer.  Testing  will  be 
conducted  by  the  Structures  Test  Branch  (WL/FIBT) , 
Mr.  John  V.  Anselmo  and  Mr.  Ray  Fisher,  Test  Engineers. 

This  report  covers  work  done  from  February  1990  through  April 
1991.  This  manuscript  was  released  by  the  author  in  May  1992  for 
publication  as  a  Wright  Laboratory  Technical  Memorandum. 
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Structural  Integrity  Branch 


This  memorandum  has  been  reviewed  and  approved. 


William  P.  Johnson,  Tech  Mgr 
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ABSTRACT 


i. 


Pressure  and  temperature  data  are  presented  for  life  cycle 
testing  of  the  B-1B  transparency.  For  each  of  three  mission  types, 
contour  plots  of  pressure  and  temperature  are  presented  for  each 
mission  segment.  Potential  flow  theory  was  used  to  calculate 
pressure.  Isentropic  flow  was  assumed  for  calculating  temperature. 
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1.  INTRODUCTION 


Aircraft  cockpit  transparencies  (canopies  or  windshields)  must 
withstand  extremes  of  temperature  and  pressure  that  result  from 
climatic  and  aerodynamic  conditions.  They  are  complex  laminates  of 
plastics  chosen  to  provide  optical  clarity  and  protection  from  bird 
impact  and  other  hazards.  They  are  thick  enough  to  generate 
significant  internal  stresses  when  the  temperature  changes  rapidly, 
such  as  when  the  aircraft  changes  speed  or  altitude. 

Life  cycle  testing  is  accomplished  by  subjecting  the 
transparency  to  a  schedule  of  temperature  and  pressure  combinations 
that  are  expected  to  produce  damage  equivalent  to  the  service  life 
of  the  transparency.  The  temperature  and  pressure  are  determined 
from  the  intended  use  of  the  aircraft  which  is  described  by  mission 
types  and  mission  segments.  Most  mission  names  are  descriptive, 
such  as  Training,  Air-To-Air  Combat,  Air-To-Ground  Combat,  etc. 
other  names,  such  as  Red  Flag,  are  less  descriptive.  Each  mission 
type  consists  of  a  number  of  mission  segments,  such  as  Takeoff, 
Cruise,  Terrain  Following,  etc.,  which  describe  the  operation  of 
the  aircraft  for  a  period  of  time  during  the  complete  mission. 
Associated  with  each  mission  segment  is  a  combination  of  aircraft 
configuration  (weight,  center  of  gravity,  wing  sweep,  flap 
position,  gear  position,  etc)  and  flight  condition  (altitude,  Mach 
number.)  This  information  is  required  to  calculate  flight  loads, 
including  pressure  and  temperature.  When  details  of  the  aircraft 
configuration  or  flight  condition  are  not  supplied,  they  must  be 
estimated  by  the  analyst. 

Selecting  mission  types  and  mission  segments  to  represent  the 
service  life  of  a  transparency  is  an  art  subject  to  much 
speculation.  It  is  a  balance  between  complexity  that  may  offer 
improved  accuracy  and  simplicity  that  will  expedite  testing.  With 
respect  to  this  technical  memorandum,  the  mission  types  and  mission 
segments  were  given.  The  objective  of  the  effort  reported  herein 
was  to  produce  compatible  temperatures  and  pressures  on  the  surface 
of  the  transparencies.  A  similar  analysis  was  accomplished  for  the 
F-lllA  aircraft  in  1987  by  Moon  [1]. 
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2 .  ANALYSIS 


s 


2.1  B-l  Fuselage  Surface  Geometry 

The  geometry  required  for  the  aerodynamic  analysis  was 
obtained  from  a  NASTRAN  model  of  the  B-1A  aircraft  forward 
fuselage.  This  was  a  computer  file  defining  the  location  of 
structural  nodes  in  three  dimensional  space.  The  NASTRAN  model  was 
originally  generated  by  ASIAC  (Aerospace  structures  Information  and 
Analysis  Center)  for  the  B-1A  aircraft  using  surface  geometry.  In 
general,  NASTRAN  models  the  structural  character  of  the  aircraft 
utilizing  geometry  interior  to  the  mold  line.  In  this  case, 
however,  ASIAC  used  the  surface  geometry,  thus  providing  the 
geometry  needed  for  an  aerodynamic  analysis.  A  graphics 
workstation  was  used  to  check  the  data.  A  few  interior  points  were 
identified  and  eliminated. 

The  surface  geometry  of  the  NASTRAN  model  was  compared  with  a 
computer  printout  for  the  B-1B  aircraft  obtained  from  the  B-l 
System  Program  Office  (SPO) .  The  most  significant  deviation  for 
the  area  of  interest  was  approximately  one  inch  located  at  a  point 
on  the  top  centerline  of  the  aircraft  near  the  aft  of  the 
transparency  (body  station  317) .  Closer  to  the  nose  the  deviation 
was  less  than  0.1  inch.  The  difference  in  geometry  between  the  B- 
1A  and  B-1B  aircraft  was  judged  to  be  insignificant  with  respect  to 
calculating  pressure  and  temperature  on  the  transparency. 


2.2  Aerodynamic  Panels 

The  surface  geometry  was  input  to  the  Configuration  Data 
Management  System  (CDMS)  [2].  CDMS  accepts  data  from  a  variety  of 
sources  and  has  a  host  of  capabilities  including  the  ability  to 
manipulate  surface  grids  while  retaining  the  surface  geometry.  The 
surface  geometry  was  input  to  CDMS  using  an  interface  program, 
TAPEIGES,  which  is  part  of  CDMS.  Utilizing  the  I3G  program  of 
CDMS,  a  surface  grid  was  generated  for  the  aerodynamic  computer 
program  QUADPAN  [3].  QUADPAN  is  a  Lockheed  Company  proprietary 
aerodynamic  panel  code. 

A  total  of  1914  aerodynamic  panels  were  used  to  modeled  one- 
half  of  the  forward  fuselage.  Body  panels  were  used  from  the  nose 
to  a  point  225  inches  aft  of  the  transparency  where  wake  panels 
extended  the  cross-section  to  a  point  1179  inches  aft  of  the 
transparency.  Of  the  1914  aerodynamic  panels,  861  were  body  panels 
located  on  the  transparency,  873  were  body  panels  not  located  on 
the  transparency  and  180  were  wake  panels.  This  paneling  was 
judged  to  provide  accurate  modeling  for  an  aerodynamic  analysis. 
A  listing  of  the  area  associated  with  each  aerodynamic  panel  in  the 
vicinity  of  the  transparency  and  the  location  of  the  aerodynamic 
center  of  each  panel  is  listed  in  appendix  D.  Appendix  D  also 
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presents  the  center  and  radius  of  curvature  for  cross-sections  in 
the  vicinity  of  the  transparency. 


2.3  Aerodynamic  Analysis 


The  CDMS  computer 
program  was  used  to 
generate  the  geometric 
input  data  for  the  QUADPAN 
aerodynamic  computer  code. 

Additional  data  was  then 
added  to  execute  QUADPAN 
for  the  combinations  of 
Mach  number  and  angles  of 
attack  listed  in  table  1. 

This  produced  an 
aerodynamic  database  for 
subsequent  use  in  an 
analysis  program  written  for  this  effort. 


Table  I 


Mach  Number 

Body  Angle- 

-of -Attack 

.30 

0 

4 

8 

.56 

0 

2 

4 

.70 

0 

1 

2  3 

.75 

0 

1 

2 

.85 

0 

1 

2 

.91 

0 

1 

2 

QUADPAN  is  a  low  order  panel  method  that  solves  the  Prandtl- 
Glauert  equation.  The  flow  field  is  assumed  to  be  derivable  from 
a  velocity  potential,  which  does  not  permit  recirculation  or 
separation  of  the  flow  field.  Shocks  are  not  considered  although 
the  Prandtl-Glauert  equation  accounts  for  Mach  number  by  geometric 
scaling  in  the  flow  direction.  The  body  is  modeled  by  a  closed 
system  of  body  panels  to  which  an  open  system  of  wake  panels  is 
attached.  Each  body  panel  consist  of  a  distributed  source  and 
doublet  singularity.  The  source  singularities  permits  the  volume 
of  the  fuselage  to  be  modeled.  Wake  panels  consisting  of  doublet 
singularities  are  attached  to  the  trailing  edge  of  the  last 
streamwise  body  panel.  The  wake  panel  singularities  are  of 
constant  strength  in  the  streamwise  direction  and  serve  to  assure 
that  the  wake  leaves  the  body  tangentially.  In  effect,  they  model 
the  starting  vortex,  and  therefore  must  extend  down-stream  far 
enough  for  the  starting  vortex  not  to  influence  the  flow  over  the 
body  panels.  Doublet  singularities  on  the  body  and  in  the  wake 
model  the  aerodynamic  lift. 

Executing  the  QUADPAN  computer  code  generates  a  dump  file  that 
contains  all  of  the  results  produced  by  the  program.  This  file  is 
in  a  format  that  is  described  in  the  QUADPAN  computer  manual.  A 
computer  program  CANOPY  was  written  to  read  the  dump  file  and 
extract  the  aerodynamic  pressure  coefficients  for  each  panel.  A 
listing  of  this  program  is  presented  in  appendix  G. 


2.4  Mission\Mission  Segment  Analysis 

Three  missions  entitled  "Dyess  Training  Mission  1",  Dyess 
Training  Mission  2”,  and  "Red  Flag  Training  Mission"  were 
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evaluated.  The  missions  segments  for  each  of  these  missions 
together  with  pressure  and  temperature  calculations  are  presented 
in  appendix  A,  B,  and  C  respectively. 

Computer  program  INTERP,  which  is  listed  in  appendix  E,  was 
written  to  calculate  pressure  and  temperature  on  the  transparency 
and  output  results  in  a  format  needed  for  the  life  cycle  test. 
INTERP  accessed  the  aerodynamic  database  generated  by  QUADPAN. 
Mach  number  data  closest  to  the  mission  segment  Mach  number  was 
selected.  Linear  interpolation  was  used  between  angle  of  attack 
tables.  An  atmospheric  database  was  incorporated  in  INTERP  that 
provided  the  absolute  pressure,  temperature,  density  and  speed  of 
sound  for  a  standard  atmosphere  as  a  function  of  altitude. 
Temperature  was  calculated  using  equation  1,  which  assumes 
adiabatic  flow. 


T  = 


Tm 

< 


P. 


\ 


/ 


(1) 


2.4.1  Output  Geometry 

Contour  plots  of  pressure  and  temperature  were  generated  for 
the  fuselage  surface  in  the  vicinity  of  the  transparency.  The 
abscissa  is  the  fuselage  station  and  the  ordinate  is  the  distance 
from  the  top  center-line  of  the  aircraft  measured  along  a  cross- 
section.  The  constant  pressure  for  each  contour  line  is  listed  in 
a  table  to  the  right  of  the  plot. 

The  ordinate  was  calculated  by  assuming  each  cross-section  of 
the  transparency  to  be  a  circular  arc.  Three  points  were  used  to 
determine  the  radius  and  center  of  curvature  of  each  arc.  Points 
closest  to  the  top  and  bottom  of  the  transparency  were  selected  as 
well  as  points  closest  to  the  middle.  The  location  of  the  center 
of  curvature  and  the  radius  of  the  circular  arc  are  given  in 
equations  2,  3  and  4.  These  equations  were  implemented  in  computer 
program  SURF  which  is  listed  in  appendix  F. 

v  _  z\ -z\  +  Y\ - Yi~2B( Z3-Zx) 

°~  2[M(Z3-Z1)*Y3-Y1] 


ZQ  =  MY0+B 


Ro  =  JlZ'-Zjt+iY'-Yj2 


(2) 

(3) 

(4) 
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Where : 


*l-*2 


M  =  - 

Z2-Zx 

y  =  y2^i 
2 

Z  =  A*  zi 
2 

B  =  H-M? 


2.4.2  Output  Pressure 

The  contour  plots  for  pressure  present  the  pressure 
differential  across  the  transparency  in  pounds  per  square  inch 
(PSI) .  A  positive  value  indicates  that  the  cockpit  pressure  is 
greater  than  the  external  pressure,  resulting  in  a  net  outward 
pressure.  The  external  pressure  was  calculated  using  equation  5, 
where  Cp  is  the  pressure  coefficient. 

P  =  p.+(|p^)cp  (5> 

The  internal  pressure  was  determined  from  the  pressurization 
schedule  for  the  B-lB  aircraft.  Cabin  pressure  was  maintained  at 
the  static  pressure  of  the  atmosphere  for  altitudes  from  sea  level 
to  8000  feet.  For  altitudes  above  8000  feet,  the  cabin  pressure 
was  maintained  at  10.916  PSI,  which  is  the  ambient  pressure  for 
8000  feet. 


2.4.3  Output  Temperature 

Contour  plots  of  temperature  present  the  temperature  of  the 
airstream  external  to  any  boundary  layer  on  surface  of  the 
transparency,  given  in  degrees  Fahrenheit.  Potential  flow  was 
assumed,  thus  neglecting  all  viscous  heating  of  the  boundary  layer. 
No  consideration  was  given  to  heat  flux  that  would  result  from 
temperature  gradients  through  the  transparency. 
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3 .  CONCLUSIONS 


1  Temperature  and  pressure  were  determined  for  the  B-1B 
transparency  for  the  purpose  of  life  cycle  testing.  The 
calculations  were  based  on  potential  flow  aerodynamics  and 
adiabatic  boundary  conditions.  It  is  believed  that  these 
assumptions  are  reasonable  for  the  range  of  Mach  numbers  up  to 
0.91. 


2  Pressure  was  presented  in  pounds  per  square  inch  differential 
pressure  across  the  transparency,  positive  for  a  net  outward  acting 
pressure. 

3  Temperature  was  presented  in  degrees  Fahrenheit  on  the  outside 
of  any  boundary- layer  on  the  surface  of  the  transparency.  Heating 
due  to  viscous  action  was  not  considered.  No  heat  transfer  was 
considered. 

4  Temperature  and  pressure  were  presented  in  the  form  of  contour 
plots  with  the  abscissa  given  in  fuselage  station  and  the  ordinate 
given  in  terms  of  distance  along  a  cross-section  from  the  top 
centerline  of  the  aircraft. 
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APPENDIX  A 


DYESS  TRAINING  MISSION  1 


DYESS  TRAINING  MISSION  #1 


LEG  WING  SWEEP 

(DEGREES) 

MACH 

ALTITUDE 
(1000  FT) 

ALPHA 

(DEG) 

TAKEOFF 

20 

.3 

2 

8 

CLIMB 

25 

.56 

10 

2 

CRUISE 

25 

.67 

20 

2 

CRUISE 

25 

.72 

20 

1.3 

DESCEND 

67.5 

.85 

9 

3 

TERRAIN  FOLLOW 

67.5 

.85 

8 

2 

CLIMB /CRUISE 

25 

.72 

16 

.7 

DESCEND 

67.5 

.85 

9 

3 

TERRAIN  FOLLOW 

67.5 

.85 

8 

1.6 

CLIMB/ CRUISE 

25 

.72 

13 

.5 

JET  PENETRATION 

25 

.72 

2 

4 

TOUCH  &  GO /LAND 

20 

.3 

2 
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TRANSPARENCY  GEOMETRY 


CURVATURE  DATA 
BODY  STATION  Y0  Z„ 


315.5750 

312.7263 
309.8799 
307.0359 

304.1943 
301.3551 
298.5182 
295.6837 
292.8515 
290.0218 

287.1944 
284.3693 
281.5467 

278.7264 
275.9085 
273.0930 
270.2798 
267.4690 
264.6605 
261.8540 
259.0495 
256.2470 
253.4463 
250.6477 
247.8509 
245.0561 
242.2633 
239.4724 
236.6834 
233.8964 
231.1113 
228.3282 
225.5470 
222.7678 
219.9906 
217.2152 
214.4418 
211.6704 
208.9009 
206.1333 


-6.7476 

-5.8132 

-4.8737 

-3.9302 

-2.9848 

-2.0393 

-1.0963 

-.1581 

.7725 

1.6922 

2.5981 

3.4861 

4.3521 

5.1916 

6.0000 

6.7722 

7.5026 

8.1863 

8.8261 

9.4257 

9.9808 

10.4868 

10.9403 

11.3366 

11.6699 

11.9373 

12.1333 

12.2551 

12.2982 

12.2604 

12.1386 

11.9323 

11.6404 

11.2621 

10.8003 

10.2559 

9.6326 

8.9346 

8.1676 

7.3373 


69.4307 

68.2033 

66.9529 

65.6796 

64.3847 

63.0693 

61.7345 

60.3820 

59.0136 

57.6314 

56.2376 

54.8349 

53.4263 

52.0152 

50.6053 

49.2006 

47.8055 

46.4243 

45.0571 

43.7024 

42.3641 

41.0462 

39.7522 

38.4862 

37.2539 

36.0586 

34.9058 

33.7986 

32.7426 

31.7406 

30.7970 

29.9136 

29.0936 

28.3391 

27.6496 

27.0264 

26.4684 

25.9740 

25.5401 

25.1637 


Ro 

37.7934 
38.4642 
39.1532 
39.8595 
40.5807 
41.3147 
42.0590 
42.8105 
43.5661 
44.3219 
45.0742 
45.8184 
46.5492 
47.2612 
47.9484 
48.6045 
49.2224 
49.7955 
50.3261 
50.8175 
51.2646 
51.6612 
52.0024 
52.2821 
52.4923 
52.6286 
52.6834 
52 . 6525 
52.5291 
52.3095 
51.9889 
51.5654 
51.0363 
50.3998 
49.6581 
48.8113 
47.8628 
46.8167 
45.6794 
44.4572 


D-l 


SURFACE  GEOMETRY  FOR  B1  TRANSPARENCY 


I  J  X 


1 

1 

315.5051 

2 

1 

315.5094 

3 

1 

315.5149 

4 

1 

315.5228 

5 

1 

315.5339 

6 

1 

315.5468 

7 

1 

315.5585 

8 

1 

315.5673 

9 

1 

315.5728 

10 

1 

315.5750 

11 

1 

315.5746 

12 

1 

315.5726 

13 

1 

315.5701 

14 

1 

315.5676 

15 

1 

315.5658 

16 

1 

315.5650 

17 

1 

315.5658 

18 

1 

315.5677 

19 

1 

315.5704 

20 

1 

315.5735 

1 

2 

312.5202 

2 

2 

312.5329 

3 

2 

312.5492 

4 

2 

312.5724 

5 

2 

312.6052 

6 

2 

312.6432 

7 

2 

312.6775 

8 

2 

312.7036 

9 

2 

312.7198 

10 

2 

312.7263 

11 

2 

312.7249 

12 

2 

312.7191 

13 

2 

312.7117 

14 

2 

312.7043 

15 

2 

312.6990 

16 

2 

312.6967 

17 

2 

312.6989 

18 

2 

312.7047 

19 

2 

312.7125 

20 

2 

312.7218 

1 

3 

309.5455 

2 

3 

309.5661 

3 

3 

309.5925 

4 

3 

309.6302 

5 

3 

309.6834 

6 

3 

309.7451 

7 

3 

309.8007 

Y  Z 

-1.8261  106.9023 

-5.4789  106.9104 

-9.1128  106.5309 

-12.7033  105.8444 

-16.2566  104.9895 

-19.7778  104.0201 

-23.2370  102.8449 

-26.5653  101.3501 

-29.6912  99.4424 

-32.4981  97.0939 

-34.9250  94.3474 

-36.9826  91.3117 

-38.6573  88.0616 

-39.9888  84.6500 

-40.9952  81.1389 

-41.7649  77.5619 

-42.4093  73.9602 

-42.9769  70.3434 

-43.5075  66.7206 

-44.0058  63.0928 

-1.8131  106.4589 

-5.4377  106.4225 

-9.0387  106.0031 

-12.5893  105.2708 

-16.0912  104.3541 

-19.5527  103.3141 

-22.9531  102.0789 

-26.2298  100.5445 

-29.3151  98.6224 

-32.0991  96.2843 

-34.5237  93.5682 

-36.5964  90.5739 

-38.3028  87.3702 

-39.6777  84.0063 

-40.7368  80.5404 

-41.5589  77.0044 

-42.2442  73.4372 

-42.8379  69.8496 

-43,3791  66.2528 

-43.8741  62.6482 

-1.7999  105.9853 

-5.3961  105.9063 

-8.9645  105.4491 

-12.4756  104.6736 

-15.9272  103.6986 

-19.3305  102.5919 

-22.6735  101.2999 


S  AREA 

1.8268  11.0480 

5.4631  11.0446 

9.0864  11.0317 

12.6927  11.0049 

16.2609  10.9697 

19.7658  10.9475 

23.2163  10.9635 

26.6333  10.9195 

30.0588  11.0410 

33.4895  10.8552 

36.9282  10.9490 

40.3659  10.7250 

43.7905  10.7513 

47.2211  10.6202 

50.6463  10.5966 

54.0855  10.5567 

57.5279  10.5406 

60.9560  10.5162 

64.3400  10.4941 

67.6531  10.4714 

1.8138  10.9418 

5.4270  10.9424 

9.0300  10.9308 

12.6194  10.8980 

16.1770  10.8477 

19.6827  10.8279 

23.1445  10.8585 

26.5786  10.8314 

30.0213  10.9534 

33.4676  10.7818 

36.9195  10.8624 

40.3689  10.6473 

43.8048  10.6611 

47.2454  10.5337 

50.6803  10.5048 

54.1287  10.4703 

57.5826  10.4618 

61.0264  10.4397 

64.4319  10.4208 

67.7746  10.4005 

1.8005  10.8384 

5.3901  10.8425 

8.9718  10.8320 

12.5425  10.7947 

16.0868  10.7339 

19.5894  10.7161 

23.0584  10.7583 
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8 

3 

309.8430 

-25.8997 

9 

3 

309.8694 

-28.9451 

10 

3 

309.8799 

-31.7061 

11 

3 

309.8776 

-34.1274 

12 

3 

309.8683 

-36.2139 

13 

3 

309.8562 

-37.9502 

14 

3 

309.8443 

-39.3665 

15 

3 

309.8355 

-40.4758 

16 

3 

309.8318 

-41.3483 

17 

3 

309.8355 

-42.0731 

18 

3 

309.8448 

-42.6922 

19 

3 

309.8576 

-43.2442 

20 

3 

309.8726 

-43.7366 

1 

4 

306.5809 

-1.7866 

2 

4 

306.6089 

-5.3543 

3 

4 

306.6449 

-8.8900 

4 

4 

306.6961 

-12.3622 

5 

4 

306.7686 

-15.7647 

6 

4 

306.8525 

-19.1112 

7 

4 

306.9280 

-22.3983 

8 

4 

306.9857 

-25.5750 

9 

4 

307.0216 

-28.5809 

10 

4 

307.0359 

-31.3188 

11 

4 

307.0327 

-33.7362 

12 

4 

307.0201 

-35.8351 

13 

4 

307.0036 

-37.5994 

14 

4 

306.9875 

-39.0551 

15 

4 

306.9755 

-40.2123 

16 

4 

306.9705 

-41.1331 

17 

4 

306.9755 

-41.8959 

18 

4 

306.9882 

-42.5399 

19 

4 

307.0056 

-43.1028 

20 

4 

307.0260 

-43.5935 

1 

5 

303.6264 

-1.7731 

2 

5 

303.6614 

-5.3122 

3 

5 

303.7063 

-8.8154 

4 

5 

303.7701 

-12.2491 

5 

5 

303.8606 

-15.6035 

6 

5 

303.9653 

-18.8949 

7 

5 

304.0596 

-22.1275 

8 

5 

304.1315 

-25.2557 

9 

5 

304.1764 

-28.2228 

10 

5 

304.1943 

-30.9375 

11 

5 

304.1903 

-33.3501 

12 

5 

304.1745 

-35.4600 

13 

5 

304.1540 

-37.2506 

14 

5 

304.1338 

-38.7434 

15 

5 

304.1189 

-39.9463 

16 

5 

304.1126 

-40.9133 

17 

5 

304.1188 

-41.7128 

18 

5 

304.1347 

-42.3810 

19 

5 

304.1564 

-42.9549 

99.7284 

26.5059 

10.7452 

97.7932 

29.9626 

10.8656 

95.4660 

33.4217 

10.7073 

92.7801 

36.8846 

10.7760 

89.8266 

40.3444 

10.5706 

86.6685 

43.7904 

10.5735 

83.3519 

47.2400 

10.4499 

79.9307 

50.6841 

10.4164 

76.4357 

54.1409 

10.3856 

72.9030 

57.6054 

10.3833 

69.3452 

61.0634 

10.3621 

65.7748 

64.4891 

10.3457 

62.1940 

67.8594 

10.3276 

105.4814 

1.7872 

10.7373 

105.3618 

5.3525 

10.7447 

104.8687 

8.9118 

10.7352 

104.0526 

12.4623 

10.6956 

103.0228 

15.9908 

10.6270 

101.8533 

19.4868 

10.6110 

100.5080 

22.9590 

10.6614 

98.9017 

26.4158 

10.6598 

96.9549 

29.8829 

10.7779 

94.6389 

33.3523 

10.6326 

91.9829 

36.8244 

10.6898 

89.0697 

40.2927 

10.4939 

85.9567 

43.7477 

10.4872 

82.6867 

47.2054 

10.3675 

79.3099 

50.6578 

10.3303 

75.8556 

54.1224 

10.3016 

72.3577 

57.5963 

10.3049 

68.8300 

61.0672 

10.2840 

65.2867 

64.5111 

10.2700 

61.7303 

67.9071 

10.2543 

104.9473 

1.7737 

10.6390 

104.7888 

5.3142 

10.6486 

104.2620 

8.8503 

10.6407 

103.4080 

12.3789 

10.5994 

102.3269 

15.8892 

10.5261 

101.0985 

19.3753 

10.5121 

99.7032 

22.8469 

10.5683 

98.0646 

26.3092 

10.5755 

96.1074 

29.7836 

10.6893 

93.8030 

33.2605 

10.5558 

91.1767 

36.7395 

10.6034 

88.3033 

40.2146 

10.4174 

85.2347 

43.6771 

10.4030 

82.0108 

47.1418 

10.2873 

78.6779 

50.6016 

10.2470 

75.2643 

54.0730 

10.2193 

71.8014 

57.5552 

10.2261 

68.3041 

61.0375 

10.2049 

64.7883 

64.4977 

10.1926 

D-3 


20 

5 

304.1819 

-43.4447 

1 

6 

300.6821 

-1.7596 

2 

6 

300.7235 

-5.2698 

3 

6 

300.7767 

-8.7406 

4 

6 

300.8524 

-12.1363 

5 

6 

300.9595 

-15.4438 

6 

6 

301.0836 

-18.6815 

7 

6 

301.1954 

-21.8610 

8 

6 

301.2806 

-24.9418 

9 

6 

301.3338 

-27.8706 

10 

6 

301.3551 

-30.5620 

11 

6 

301.3503 

-32.9691 

12 

6 

301.3316 

-35.0887 

13 

6 

301.3073 

-36.9036 

14 

6 

301.2833 

-38.4315 

15 

6 

301.2657 

-39.6777 

16 

6 

301.2582 

-40.6889 

17 

6 

301.2656 

-41.5237 

18 

6 

301.2843 

-42.2154 

19 

6 

301.3101 

-42.8005 

20 

6 

301.3403 

-43.2902 

1 

7 

297.7479 

-1.7459 

2 

7 

297.7953 

-5.2271 

3 

7 

297.8561 

-8.6655 

4 

7 

297.9427 

-12.0238 

5 

7 

298.0653 

-15.2855 

6 

7 

298.2073 

-18.4711 

7 

7 

298.3353 

-21.5988 

8 

7 

298.4330 

-24.6333 

9 

7 

298.4938 

-27.5243 

10 

7 

298.5182 

-30.1924 

11 

7 

298.5128 

-32.5931 

12 

7 

298.4913 

-34.7211 

13 

7 

298.4634 

-36.5586 

14 

7 

298.4360 

-38.1194 

15 

7 

298.4158 

-39.4067 

16 

7 

298.4072 

-40.4599 

17 

7 

298.4157 

-41.3285 

18 

7 

298.4372 

-42.0432 

19 

7 

298.4666 

-42.6396 

20 

7 

298.5013 

-43.1301 

1 

8 

294.8239 

-1.7321 

2 

8 

294.8768 

-5.1841 

3 

8 

294.9445 

-8.5903 

4 

8 

295.0411 

-11.9116 

5 

8 

295.1780 

-15.1286 

6 

8 

295.3365 

-18.2636 

7 

8 

295.4794 

-21.3410 

8 

8 

295.5885 

-24.3302 

9 

8 

295.6565 

-27.1840 

10 

8 

295.6837 

-29.8286 

11 

8 

295.6776 

-32.2222 

61.2569 

67.9174 

10.1792 

104.3830 

1.7602 

10.5426 

104.1876 

5.2753 

10.5549 

103.6290 

8.7872 

10.5477 

102.7396 

12.2926 

10.5065 

101.6108 

15.7827 

10.4312 

100.3275 

19.2557 

10.4185 

98.8855 

22.7229 

10.4778 

97.2169 

26.1872 

10.4922 

95.2506 

29.6655 

10.6011 

92.9584 

33.1468 

10.4789 

90.3615 

36.6305 

10.5172 

87.5273 

40.1106 

10.3407 

84.5026 

43.5791 

10.3198 

81.3241 

47.0496 

10.2081 

78.0349 

50.5157 

10.1657 

74.6617 

53.9930 

10.1382 

71.2339 

57.4821 

10.1476 

67.7675 

60.9740 

10.1253 

64.2797 

64.4486 

10.1147 

60.7740 

67.8894 

10.1030 

103.7884 

1.7465 

10.4487 

103.5580 

5.2357 

10.4618 

102.9697 

8.7226 

10.4563 

102.0475 

12.2037 

10.4159 

100.8745 

15.6714 

10.3402 

99.5403 

19.1286 

10.3293 

98.0549 

22.5879 

10.3907 

96.3586 

26.0505 

10.4095 

94.3847 

29.5292 

10.5118 

92.1050 

33.0123 

10.4004 

89.5372 

36.4984 

10.4311 

86.7417 

39.9815 

10.2650 

83.7602 

43.4544 

10.2384 

80.6267 

46.9291 

10.1303 

77.3806 

50.4005 

10.0864 

74.0478 

53.8824 

10.0577 

70.6552 

57.3770 

10.0684 

67.2202 

60.8767 

10.0444 

63.7610 

64.3633 

10.0351 

60.2815 

67.8229 

10.0256 

103.1635 

1.7326 

10.3559 

102.9000 

5.1955 

10.3707 

102.2840 

8.6567 

10.3665 

101.3317 

12.1122 

10.3274 

100.1180 

15.5560 

10.2535 

98.7367 

18.9945 

10.2441 

97.2114 

22.4426 

10.3052 

95.4899 

9000 

10.3272 

93.5096 

29.3759 

10.4224 

91.2429 

32.8578 

10.3208 

88.7040 

36.3438 

10.3444 

D-4 


12 

8 

295.6536 

-34.3573 

13 

8 

295.6225 

-36.2154 

14 

8 

295.5919 

-37.8070 

15 

8 

295.5693 

-39.1331 

16 

8 

295.5597 

-40.2263 

17 

8 

295.5692 

-41.1273 

18 

8 

295.5932 

-41.8644 

19 

8 

295.6261 

-42.4723 

20 

8 

295.6648 

-42.9643 

1 

9 

291.9099 

-1.7181 

2 

9 

291.9678 

-5.1408 

3 

9 

292.0420 

-8.5148 

4 

9 

292.1478 

-11.7997 

5 

9 

292.2976 

-14.9731 

6 

9 

292.4713 

-18.0590 

7 

9 

292.6278 

-21.0876 

8 

9 

292.7473 

-24.0325 

9 

9 

292.8218 

-26.8497 

10 

9 

292.8515 

-29.4708 

11 

9 

292.8449 

-31.8563 

12 

9 

292.8186 

-33.9972 

13 

9 

292.7845 

-35.8741 

14 

9 

292.7509 

-37.4944 

15 

9 

292.7262 

-38.8571 

16 

9 

292.7157 

-39.9881 

17 

9 

292.7261 

-40.9202 

18 

9 

292.7523 

-41.6790 

19 

9 

292.7885 

-42.2985 

20 

9 

292.8309 

-42.7928 

1 

10 

289.0061 

-1.7040 

2 

10 

289.0685 

-5.0972 

3 

10 

289.1486 

-8.4392 

4 

10 

289.2626 

-11.6881 

5 

10 

289.4241 

-14.8190 

6 

10 

289.6114 

-17.8573 

7 

10 

289.7803 

-20.8385 

8 

10 

289.9092 

-23.7401 

9 

10 

289.9897 

-26.5213 

10 

10 

290.0218 

-29.1188 

11 

10 

290.0146 

-31.4956 

12 

10 

289.9862 

-33.6408 

13 

10 

289.9494 

-35.5348 

14 

10 

289.913 2 

-37.1816 

15 

10 

289.8865 

-38.5784 

16 

10 

289.8752 

-39.7453 

17 

10 

289.8864 

-40.7070 

18 

10 

289.9148 

-41.4869 

19 

10 

289.9537 

-42.1181 

20 

10 

289.9995 

-42.6156 

1 

11 

286.1125 

-1.6898 

2 

11 

286.1789 

-5.0534 

3 

11 

286.2641 

-8.3633 

85.9465 

39.8280 

10.1879 

83.0076 

43.3035 

10.1575 

79.9185 

46.7811 

10.0535 

76.7152 

50.2563 

10.0083 

73.4225 

53.7415 

9.9781 

70.0654 

57.2400 

9.9895 

66.6621 

60.7455 

9.9630 

63.2320 

64.2419 

9.9545 

59.7794 

67.7174 

9.9466 

102.5084 

1.7186 

10.2650 

102.2137 

5.1547 

10.2807 

101.5720 

8.5895 

10.2771 

100.5922 

12.0185 

10.2404 

99.3414 

15.4366 

10.1698 

97.9169 

18.8541 

10.1620 

96.3550 

22.2880 

10.2213 

94.6106 

25.7367 

10.2452 

92.6253 

29.2066 

10.3330 

90.3719 

32.6846 

10.2402 

87.8617 

36.1679 

10.2575 

85.1418 

39.6510 

10.1113 

82.2449 

43.1272 

10.0777 

79.1995 

46.6061 

9.9772 

76.0387 

50.0837 

9.9319 

‘  72.7860 

53.5708 

9.8987 

69.4645 

57.0715 

9.9098 

66.0933 

60.5807 

9.8809 

62.6928 

64.0844 

9.8727 

59.2677 

67.5729 

9.8661 

101.8230 

1.7045 

10.1756 

101.4990 

5.1134 

10.1916 

100.8336 

8.5211 

10.1892 

99.8290 

11.9227 

10.1549 

98.5446 

15.3136 

10.0885 

97.0809 

18.7079 

10.0825 

95.4856 

22.1247 

10.1389 

93.7208 

25.5615 

10.1633 

91.7318 

29.0223 

10.2430 

89.4922 

32.4934 

10.1589 

87.0104 

35.9717 

10.1701 

84.3275 

39.4515 

10.0344 

81.4719 

42.9266 

9.9980 

78.4698 

46.4049 

9.9015 

75.3511 

49.8831 

9.8564 

72.1382 

53.3706 

9.8203 

68.8525 

56.8716 

9.8302 

65.5138 

60.3823 

9.7979 

62.1434 

63.8906 

9.7899 

58.7465 

67.3891 

9.7849 

101.1074 

1.6901 

10.0871 

100.7560 

5.0715 

10.1034 

100.0690 

8.4515 

10.1016 

D-5 


4 

11 

286.3855 

-11.5768 

5 

11 

286.5575 

-14.6663 

6 

11 

286.7570 

-17.6586 

7 

11 

286.9370 

-20.5937 

8 

11 

287.0744 

-23.4532 

9 

11 

287.1601 

-26.1989 

10 

11 

287.1944 

-28.7726 

11 

11 

287.1867 

-31.1399 

12 

11 

287.1565 

-33.2881 

13 

11 

287.1173 

-35.1974 

14 

11 

287.0786 

-36.8686 

15 

11 

287.0502 

-38.2973 

16 

11 

287.0381 

-39.4979 

17 

11 

287.0501 

-40.4878 

18 

11 

287.0803 

-41.2882 

19 

11 

287.1218 

-41.9313 

20 

11 

287.1706 

-42.4327 

1 

12 

283.2290 

-1.6754 

2 

12 

283.2990 

-5.0092 

3 

12 

283.3887 

-8.2872 

4 

12 

283.5166 

-11.4658 

5 

12 

283.6978 

-14.5151 

6 

12 

283.9081 

-17.4628 

7 

12 

284.0979 

-20.3533 

8 

12 

284.2428 

-23.1717 

9 

12 

284.3332 

-25.8824 

10 

12 

284.3693 

-28.4323 

11 

12 

284.3613 

-30.7892 

12 

12 

284.3294 

-32.9391 

13 

12 

284.2880 

-34.8618 

14 

12 

284.2473 

-36.5554 

15 

12 

284.2172 

-38.0137 

16 

12 

284.2045 

-39.2458 

17 

12 

284.2171 

-40.2625 

18 

12 

284.2490 

-41.0829 

19 

12 

284.2928 

-41.7380 

20 

12 

284.3442 

-42.2441 

1 

13 

280.3556 

-1.6609 

2 

13 

280.4286 

-4.9647 

3 

13 

280.5222 

-8.2109 

4 

13 

280.6558 

-11.3551 

5 

13 

280.8451 

-14.3652 

6 

13 

281.0648 

-17.2699 

7 

13 

281.2630 

-20.1173 

8 

13 

281.4144 

-22.8955 

9 

13 

281.5089 

-25.5719 

10 

13 

281.5467 

-28.0979 

11 

13 

281.5383 

-30.4437 

12 

13 

281.5049 

-32.5939 

13 

13 

281.4617 

-34.5282 

14 

13 

281.4191 

-36.2420 

15 

13 

281.3877 

-37.7276 

99.0421 

11.8249 

10.0700 

97.7276 

15.1875 

10.0091 

96.2286 

18.5565 

10.0048 

94.6034 

21.9535 

10.0575 

92.8205 

25.3755 

10.0813 

90.8291 

28.8242 

10.1522 

88.6037 

32.2855 

10.0759 

86.1501 

35.7564 

10.0828 

83.5037 

39.2305 

9.9570 

80.6888 

42.7025 

9.9191 

77.7293 

46.1784 

9.8262 

74.6523 

49.6556 

9.7814 

71.4791 

53.1417 

9.7414 

68.2294 

56.6410 

9.7498 

64.9236 

60.1509 

9.7138 

61.5838 

63.6611 

9.7059 

58.2156 

67.1666 

9.7019 

100.3616 

1.6759 

10.0001 

99.9846 

5.0290 

10.0159 

99.2780 

8.3809 

10.0146 

98.2315 

11.7253 

9.9863 

96.8904 

15.0585 

9.9311 

95.3601 

18.4005 

9.9281 

93.7083 

21.7752 

9.9764 

91.9096 

25.1797 

9.9989 

89.9172 

28.6132 

10.0609 

87.7064 

32.0622 

9.9922 

85.2808 

35.5229 

9.9949 

82.6702 

38.9893 

9.8784 

79.8955 

42.4558 

9.8399 

76.9780 

45.9276 

9.7513 

73.9424 

49.4018 

9.7069 

70.8087 

52.8847 

9.6628 

67.5952 

56.3802 

9.6691 

64.3226 

59.8870 

9.6296 

61.0140 

63.3963 

9.6208 

57.6752 

66.9053 

9.6175 

99.5855 

1.6613 

9.9137 

99.1848 

4.9861 

9.9290 

98.4607 

8.3092 

9.9281 

97.3971 

11.6241 

9.9025 

96.0331 

14.9269 

9.8535 

94.4753 

18.2404 

9.8520 

92.8002 

21.5907 

9.8956 

90.9882 

24.9748 

9.9160 

88.9961 

28.3905 

9.9692 

86.8003 

31.8244 

9.9072 

84.4024 

35.2727 

9.9063 

81.8272 

38.7290 

9.7997 

79.0920 

42.1880 

9.7607 

76.2160 

45.6534 

9.6758 

73.2213 

49.1228 

9.6327 

D-6 


16 

13 

281.3744 

17 

13 

281.3876 

18 

13 

281.4209 

19 

13 

281.4667 

20 

13 

281.5204 

1 

14 

277.4923 

2 

14 

277.5679 

3 

14 

277.6649 

4 

14 

277.8031 

5 

14 

277.9992 

6 

14 

278.2268 

7 

14 

278.4323 

8 

14 

278.5892 

9 

14 

278.6873 

10 

14 

278.7264 

11 

14 

278.7177 

12 

14 

278.6831 

13 

14 

278.6382 

14 

14 

278.5941 

15 

14 

278.5616 

16 

14 

278.5478 

17 

14 

278.5614 

18 

14 

278.5959 

19 

14 

278.6434 

20 

14 

278.6992 

1 

15 

274.6393 

2 

15 

274.7169 

3 

15 

274.8165 

4 

15 

274.9586 

5 

15 

275.1602 

6 

15 

275.3943 

7 

15 

275.6057 

8 

15 

275.7673 

9 

15 

275.8681 

10 

15 

275.9085 

11 

15 

275.8995 

12 

15 

275.8639 

13 

15 

275.8177 

14 

15 

275.7722 

15 

15 

275.7388 

16 

15 

275.7246 

17 

15 

275.7386 

18 

15 

275.7742 

19 

15 

275.8231 

20 

15 

275.8805 

1 

16 

271.7963 

2 

16 

271.8755 

3 

16 

271.9772 

4 

16 

272.1223 

5 

16 

272.3282 

6 

16 

272.5673 

7 

16 

272.7834 

-38.9893  70.1271 

-40.0313  66.9498 

-40.8709  63.7109 

-41.5382  60.4341 

-42.0499  57.1251 

-1.6463  98.7791 

-4.9199  98.3567 

-8.1344  97.6171 

-11.2446  96.5391 

-14.2168  95.1556 

-17.0800  93.5743 

-19.8856  91.8793 

-22.6248  90.0563 

-25.2673  88.0658 

-27.7693  85.8855 

-30.1032  83.5151 

-32.2524  80.9746 

-34.1964  78.2783 

-35.9282  75.4432 

-37.4389  72.4891 

-38.7280  69.4342 

-39.7941  66.2933 

-40.6523  63.0885 

-41.3319  59.8438 

-41.8500  56.5655 

-1.6315  97.9425 

-4.8748  97.5002 

-8.0578  96.7472 

-11.1345  95.6573 

-14.0698  94.2579 

-16.8930  92.6570 

-19.6583  90.9455 

-22.3595  89.1139 

-24.9687  87.1264 

-27.4466  84.9619 

-29.7678  82.6187 

-31.9146  80.1125 

-33.8666  77.4545 

-35.6143  74.6597 

-37.1478  71.7457 

-38.4622  68.7299 

-39.5509  65.6256 

-40.4271  62.4553 

-41.1191  59.2434 

-41.6444  55.9963 

-1.6166  97.0757 

-4.8295  96.6154 

-7.9809  95.8510 

-11.0246  94.7519 

-13.9242  93.3400 

-16.7090  91.7234 

-19.4353  89.9987 


52.6005  9.5843 

56.0901  9.5877 

59.5912  9.5441 

63.0965  9.5348 

66.6060  9.5321 

1.6467  9.8277 

4.9426  9.8424 

8.2366  9.8418 

11.5213  9.8188 

14.7931  9.7764 

18.0768  9.7766 

21.4006  9.8148 

24.7620  9.8326 

28.1572  9.8772 

31.5735  9.8216 

35.0068  9.8172 

38.4507  9.7201 

41.9000  9.6815 

45.3570  9.6003 

48.8194  9.5582 

52.2898  9.5055 

55.7715  9.5059 

59.2642  9.4581 

62.7630  9.4476 

66.2691  9.4451 

1.6318  9.7426 

4.8986  9.7560 

8.1632  9.7552 

11.4173  9.7350 

14.6574  9.6989 

17.9100  9.7003 

21.2057  9,7333 

24.5421  9.7480 

27.9144  9.7839 

31.3107  9.7343 

34.7268  9.7274 

38.1559  9.6391 

41.5932  9.6014 

45.0396  9.5238 

48.4932  9.4835 

51.9542  9.4263 

55.4255  9.4232 

58.9074  9.3708 

62.3964  9.3592 

65.8959  9.3571 

1.6169  9.6578 

4.8542  9.6699 

8.0888  9.6685 

11.3119  9.6507 

14.5199  9.6205 

17.7407  9.6236 

21.0066  9.6513 


8 

16 

272.9485 

9 

16 

273.0517 

10 

16 

273.0930 

11 

16 

273.0838 

12 

16 

273.0474 

13 

16 

273.0002 

14 

16 

272.9537 

15 

16 

272.9194 

16 

16 

272.9050 

17 

16 

272.9193 

18 

16 

272.9557 

19 

16 

273.0056 

20 

16 

273.0644 

1 

17 

268.9634 

2 

17 

269.0437 

3 

17 

269.1469 

4 

17 

269.2940 

5 

17 

269.5030 

6 

17 

269.7458 

7 

17 

269.9652 

8 

17 

270.1330 

9 

17 

270.2379 

10 

17 

270.2798 

11 

17 

270.2705 

12 

17 

270.2335 

13 

17 

270.1855 

14 

17 

270.1382 

15 

17 

270.1034 

16 

17 

270.0887 

17 

17 

270.1033 

18 

17 

270.1403 

19 

17 

270.1910 

20 

17 

270.2507 

1 

18 

266.1406 

2 

18 

266.2215 

3 

18 

266.3255 

4 

18 

266.4739 

5 

18 

266.6847 

6 

18 

266.9297 

7 

18 

267.1513 

8 

18 

267.3207 

9 

18 

267.4266 

10 

18 

267.4690 

11 

18 

267.4596 

12 

18 

267.4222 

13 

18 

267.3737 

14 

18 

267.3260 

15 

18 

267.2908 

16 

18 

267.2760 

17 

18 

267.2907 

18 

18 

267.3280 

19 

18 

267.3793 

-22.0995  88.1609 

-24.6761  86.1777 

-27.1298  84.0295 

-29.4375  81.7133 

-31.5806  79.2408 

-33.5386  76.6204 

-35.3002  73.8654 

-36.8541  70.9912 

-38.1918  68.0143 

-39.3016  64.9469 

-40.1952  61.8115 

-40.8998  58.6328 

-41.4331  55.4174 

-1.6015  96.1786 

-4.7838  95.7022 

-7.9037  94.9284 

-10.9151  93.8227 

-13.7800  92.4019 

-16.5278  90.7736 

-19.2167  89.0390 

-21.8449  87.1975 

-24.3894  85.2198 

-26.8188  83.0884 

-29.1122  80.7989 

-31.2503  78.3595 

-33.2126  75.7761 

-34.9859  73.0603 

-36.5579  70.2256 

-37.9167  67.2875 

-39.0463  64.2570 

-39.9567  61.1568 

-40.6740  58.0120 

-41.2161  54.8291 

-1.5864  95.2515 

-4.7378  94.7609 

-7.8264  93.9797 

-10.8059  92.8700 

-13.6372  91.4438 

-16.3496  89.8076 

-19.0023  88.0665 

-21.5958  86.2234 

-24.1087  84.2528 

-26.5138  82.1385 

-28.7920  79.8754 

-30.9237  77.4687 

-32.8884  74.9217 

-34.6713  72.2445 

-36.2592  69.4488 

-37.6371  66.5494 

-38.7851  63.5560 

-39.7116  60.4916 

-40.4418  57.3810 


24.3160  9.6632 

27.6634  9.6904 

31.0373  9.6460 

34.4338  9.6371 

37.8460  9.5575 

41.2691  9.5210 

44.7027  9.4467 

48.1450  9.4085 

51.5946  9.3465 

55.0531  9.3396 

58.5216  9.2834 

61.9981  9.2703 

65.4877  9.2680 

1.6017  9.5732 

4.8093  9.5840 

8.0136  9.5818 

11.2056  9.5659 

14.3809  9.5410 

17.5691  9.5454 

20.8041  9.5681 

24.0846  9.5771 

27.4051  9.5964 

30.7544  9.5568 

34.1292  9.5463 

37.5222  9.4749 

40.9289  9.4398 

44.3476  9.3691 

47.7762  9.3325 

51.2121  9.2664 

54.6557  9.2556 

58.1084  9.1949 

61.5692  9.1800 

65.0455  9.1775 

1.5867  9.4886 

4.7639  9.4976 

7.9377  9.4947 

11.0983  9.4801 

14.2408  9.4601 

17.3958  9.4662 

20.5987  9.4842 

23.8492  9.4896 

27.1406  9.5014 

30.4635  9.4664 

33.8145  9.4544 

37.1861  9.3907 

40.5741  9.3576 

43.9756  9.2902 

47.3886  9.2555 

50.8085  9.1851 

54.2348  9.1706 

57.6691  9.1055 

61.1118  9.0893 


20 

18 

267.4396 

-40.9935 

54.2311 

64.5715 

1 

19 

263.3268 

-1.5711 

94.2974 

1.5714 

2 

19 

263.4081 

-4.6916 

93.7941 

4.7182 

3 

19 

263.5123 

-7.7489 

93.0071 

7.8612 

4 

19 

263.6612 

-10.6969 

91.8956 

10.9902 

5 

19 

263.8727 

-13.4957 

90.4670 

14.0998 

6 

19 

264.1188 

-16.1742 

88.8264 

17.2216 

7 

19 

264.3412 

-18.7922 

87.0818 

20.3917 

8 

19 

264.5115 

-21.3517 

85.2394 

23.6105 

9 

19 

264.6179 

-23.8336 

83.2769 

26.8716 

10 

19 

264.6605 

-26.2143 

81.1802 

30.1663 

11 

19 

264.6510 

-28.4767 

78.9434 

33.4916 

12 

19 

264.6134 

-30.6007 

76.5687 

36.8400 

13 

19 

264.5647 

-32.5661 

74.0576 

40.2073 

14 

19 

264.5168 

-34.3565 

71.4185 

43.5895 

15 

19 

264.4814 

-35.9581 

68.6615 

46.9847 

16 

19 

264.4665 

-37.3531 

65.8006 

50.3865 

17 

19 

264.4814 

-38.5181 

62.8445 

53.7935 

18 

19 

264.5188 

-39.4602 

59.8161 

57.2073 

19 

19 

264.5703 

-40.2034 

56.7403 

60.6293 

20 

19 

264.6309 

-40.7654 

53.6240 

64.0696 

1 

20 

260.5212 

-1.5557 

93.3192 

1.5560 

2 

20 

260.6023 

-4.6451 

92.8045 

4.6721 

3 

20 

260.7065 

-7.6713 

92.0133 

7.7843 

4 

20 

260.8552 

-10.5882 

90.9019 

10.8815 

5 

20 

261.0665 

-13.3554 

89.4736 

13.9583 

6 

20 

261.3124 

-16.0012 

87.8316 

17.0465 

7 

20 

261.5347 

-18.5857 

86.0861 

20.1833 

8 

20 

261.7050 

-21.1124 

84.2464 

23.3697 

9 

20 

261.8115 

-23.5637 

82.2932 

26.5990 

10 

20 

261.8540 

-25.9198 

80.2142 

29.8640 

11 

20 

261.8445 

-28.1657 

78.0035 

33.1620 

12 

20 

261.8069 

-30.2809 

75.6604 

36.4856 

13 

20 

261.7582 

-32.2454 

73.1846 

39.8303 

14 

20 

261.7103 

-34.0415 

70.5832 

43.1913 

15 

20 

261.6750 

-35.6549 

67.8646 

46.5668 

16 

20 

261.6600 

-37.0652 

65.0421 

49.9488 

17 

20 

261.6748 

-38.2460 

62.1234 

53.3345 

18 

20 

261.7123 

-39.2031 

59.1313 

56.7260 

19 

20 

261.7639 

-39.9593 

56.0907 

60.1252 

20 

20 

261.8244 

-40.5324 

53.0085 

63.5436 

1 

21 

257.7234 

-1.5402 

92.3173 

1.5404 

2 

21 

257.8041 

-4.5985 

91.7928 

4.6258 

3 

21 

257.9077 

-7.5935 

90.9988 

7.7068 

4 

21 

258.0556 

-10.4796 

89.8895 

10.7721 

5 

21 

258.2658 

-13.2162 

88.4641 

13.8162 

6 

21 

258.5103 

-15.8306 

86.8238 

16.8708 

7 

21 

258.7317 

-18.3828 

85.0799 

19.9739 

8 

21 

258.9011 

-20.8774 

83.2447 

23.1269 

9 

21 

259.0071 

-23.2987 

81.3019 

26.3234 

10 

21 

259.0495 

-25.6302 

79.2411 

29.5574 

11 

21 

259.0401 

-27.8589 

77.0562 

32.8266 

D-9 


9.0858 

9.4041 

9.4112 

9.4076 

9.3939 

9.3780 

9.3853 

9.3992 

9.4021 

9.4061 

9.3752 

9.3623 

9.3062 

9.2749 

9.2108 

9.1780 

9.1035 

9.0854 

9.0160 

8.9978 

8.9937 

9.3193 

9.3248 

9.3200 

9.3070 

9.2952 

9.3039 

9.3131 

9.3134 

9.3107 

9.2838 

9.2704 

9.2209 

9.1917 

9.1305 

9.1001 

9.0220 

9.0001 

8.9269 

8.9065 

8.9013 

9.2342 

9.2379 

9.2316 

9.2194 

9.2106 

9.2207 

9.2264 

9.2240 

9.2150 

9.1910 

9.1766 


12 

21 

259.0026 

-29.9642 

13 

21 

258.9541 

-31.9262 

14 

21 

258.9064 

-33.7264 

15 

21 

258.8712 

-35.3495 

16 

21 

258.8564 

-36.7734 

17 

21 

258.8711 

-37.9689 

18 

21 

258.9084 

-38.9406 

19 

21 

258.9597 

-39.7100 

20 

21 

259.0201 

-40.2948 

1 

22 

254.9337 

-1.5245 

2 

22 

255.0134 

-4.5515 

3 

22 

255.1160 

-7.5155 

4 

22 

255.2624 

-10.3714 

5 

22 

255.4705 

-13.0781 

6 

22 

255.7128 

-15.6623 

7 

22 

255.9320 

-18.1834 

8 

22 

256.0999 

-20.6468 

9 

22 

256.2049 

-23.0386 

10 

22 

256.2470 

-25.3455 

11 

22 

256.2376 

-27.5562 

12 

22 

256.2005 

-29.6504 

13 

22 

256.1524 

-31.6086 

14 

22 

256.1051 

-33.4110 

15 

22 

256.0703 

-35.0422 

16 

22 

256.0556 

-36.4780 

17 

22 

256.0702 

-37.6870 

18 

22 

256.1071 

-38.6726 

19 

22 

256.1580 

-39.4553 

20 

22 

256.2177 

-40.0525 

1 

23 

252.1518 

-1.5088 

2 

23 

252.2303 

-4.5043 

3 

23 

252.3314 

-7.4373 

4 

23 

252.4757 

-10.2633 

5 

23 

252.6808 

-12.9412 

6 

23 

252.9196 

-15.4964 

7 

23 

253.1357 

-17.9875 

8 

23 

253.3013 

-20.4206 

9 

23 

253.4049 

-22.7834 

10 

23 

253.4463 

-25.0655 

11 

23 

253.4371 

-27.2577 

12 

23 

253.4005 

-29.3398 

13 

23 

253.3531 

-31.2926 

14 

23 

253.3065 

-33.0955 

15 

23 

253.2721 

-34.7328 

16 

23 

253.2576 

-36.1788 

17 

23 

253.2720 

-37.4002 

18 

23 

253.3084 

-38.3993 

19 

23 

253.3585 

-39.1954 

20 

23 

253.4175 

-39.8055 

1 

24 

249.3778 

-1.4930 

2 

24 

249.4549 

-4.4569 

3 

24 

249.5539 

-7.3590 

74.7443 

36.1236 

9.1340 

72.3034 

39.4439 

9.1073 

69.7391 

42.7822 

9.0495 

67.0586 

46.1360 

9.0207 

64.2744 

49.4963 

8.9392 

61.3932 

52.8591 

8.9137 

58.4378 

56.2263 

8.8366 

55.4328 

59.6008 

8.8140 

52.3852 

62.9949 

8.8078 

91.2917 

1.5247 

9.1488 

90.7587 

4.5790 

9.1506 

89.9636 

7.6288 

9.1431 

88.8584 

10.6624 

9.1310 

87.4385 

13.6737 

9.1248 

85.8028 

16.6948 

9.1358 

84.0633 

19.7638 

9.1382 

82.2345 

22.8827 

9.1331 

80.3030 

26.0455 

9.1186 

78.2608 

29.2472 

9.0976 

76.1015 

32.4861 

9.0830 

73.8204 

35.7548 

9.0463 

71.4138 

39.0494 

9.0219 

68.8862 

42.3629 

8.9666 

66.2435 

45.6933 

8.9402 

63.4975 

49.0302 

8.8560 

60.6540 

52.3683 

8.8267 

57.7355 

55.7096 

8.7466 

54.7666 

59.0573 

8.7217 

51.7541 

62.4250 

8.7137 

90.2425 

1.5089 

9.0628 

89.7025 

4.5319 

9.0628 

88.9077 

7.5503 

9.0538 

87.8085 

10.5521 

9.0411 

86.3969 

13.5309 

9.0368 

84.7689 

16.5186 

9.0491 

83.0362 

19.5533 

9.0487 

81.2157 

22.6376 

9.0412 

79.2967 

25.7659 

9.0218 

77.2733 

28.9341 

9.0031 

75.1395 

32.1415 

8.9885 

72.8886 

35.3807 

8.9575 

70.5159 

38.6475 

8.9351 

68.0246 

41.9347 

8.8828 

65.4194 

45.2399 

8.8585 

62.7113 

48.5516 

8.7719 

59.9058 

51.8633 

8.7387 

57.0246 

55.1770 

8.6555 

54.0920 

58.4965 

8.6287 

51.1151 

61.8357 

8.6193 

89.1697 

1.4932 

8.9764 

88.6240 

4.4846 

8.9743 

87.8311 

7.4714 

8.9640 

D-10 


4 

24 

249.6954 

-10.1556 

5 

24 

249.8966 

-12.8054 

6 

24 

250.1308 

-15.3328 

7 

24 

250.3428 

-17.7952 

8 

24 

250.5053 

-20.1988 

9 

24 

250.6070 

-22.5330 

10 

24 

250.6477 

-24.7903 

11 

24 

250.6386 

-26.9634 

12 

24 

250.6026 

-29.0322 

13 

24 

250.5561 

-30.9781 

14 

24 

250.5104 

-32.7798 

15 

24 

250.4766 

-34.4213 

16 

24 

250.4624 

-35.8759 

17 

24 

250.4765 

-37.1085 

18 

24 

250.5123 

-38.1205 

19 

24 

250.5615 

-38.9301 

20 

24 

250.6194 

-39.5540 

1 

25 

246.6118 

-1.4771 

2 

25 

246.6869 

-4.4093 

3 

25 

246.7835 

-7.2805 

4 

25 

246.9215 

-10.0480 

5 

25 

247.1178 

-12.6707 

6 

25 

247.3464 

-15.1716 

7 

25 

247.5533 

-17.6063 

8 

25 

247.7120 

-19.9814 

9 

25 

247.8112 

-22.2875 

10 

25 

247.8509 

-24.5200 

11 

25 

247.8420 

-26.6733 

12 

25 

247.8070 

-28.7277 

13 

25 

247.7616 

-30.6651 

14 

25 

247.7169 

-32.4639 

15 

25 

247.6839 

-34.1078 

16 

25 

247.6700 

-35.5692 

17 

25 

247.6838 

-36.8119 

18 

25 

247.7188 

-37.8363 

19 

25 

247.7668 

-38.6596 

20 

25 

247.8233 

-39.2977 

1 

26 

243.8536 

-1.4611 

2 

26 

243.9265 

-4.3615 

3 

26 

244.0202 

-7.2019 

4 

26 

244.1541 

-9.9407 

5 

26 

244.3446 

-12.5372 

6 

26 

244.5664 

-15.0128 

7 

26 

244.7673 

-17.4210 

8 

26 

244.9212 

-19.7684 

9 

26 

245.0176 

-22.0469 

10 

26 

245.0561 

-24.2544 

11 

26 

245.0475 

-26.3873 

12 

26 

245.0135 

-28.4263 

13 

26 

244.9694 

-30.3538 

14 

26 

244.9259 

-32.1479 

15 

26 

244.8940 

-33.7923 

86.7398 

10.4416 

8.9505 

85.3392 

13.3880 

8.9475 

83.7218 

16.3424 

8.9607 

81.9987 

19.3429 

8.9581 

80.1883 

22.3921 

8.9485 

78.2828 

25.4850 

8.9242 

76.2786 

28.6189 

8.9078 

74.1700 

31.7936 

8.8928 

71.9491 

35.0016 

8.8672 

69.6097 

38.2393 

8.8469 

67.1542 

41.4984 

8.7974 

64.5862 

44.7766 

8.7751 

61.9160 

48.0615 

8.6862 

59.1484 

51.3453 

8.6506 

56.3048 

54.6299 

8.5642 

53.4092 

57.9193 

8.5351 

50.4682 

61.2285 

8.5239 

88.0732 

1.4773 

8.8895 

87.5233 

4.4370 

8.8854 

86.7337 

7.3921 

8.8736 

85.6523 

10.3306 

8.8589 

84.2654 

13.2451 

8.8561 

82.6617 

16.1666 

8.8698 

80.9507 

19.1327 

8.8653 

79.1523 

22.1467 

8.8543 

77.2614 

25.2036 

8.8260 

75.2767 

28.3023 

8.8119 

73.1932 

31.4432 

8.7966 

71.0016 

34.6188 

8.7755 

68.6952 

37.8256 

8.7574 

66.2750 

41.0551 

8.7105 

63.7439 

44.3047 

8.6902 

61.1116 

47.5609 

8.5997 

58.3820 

50.8154 

8.5610 

55.5763 

54.0695 

8.4722 

52.7179 

57.3276 

8.4411 

49.8135 

60.6048 

8.4286 

86.9531 

1.4612 

8.8018 

86.4004 

4.3891 

8.7960 

85.6157 

7.3125 

8.7824 

84.5461 

10.2194 

8.7661 

83.1755 

13.1022 

8.7630 

81.5885 

15.9912 

8.7770 

79.8923 

18.9231 

8.7713 

78.1078 

21.9015 

8.7590 

76.2324 

24.9222 

8.7274 

74.2676 

27.9848 

8.7146 

72.2089 

31.0909 

8.6995 

70.0464 

34.2329 

8.6822 

67.7723 

37.4075 

8.6663 

65.3870 

40.6058 

8.6219 

62.8925 

43.8250 

8.6031 
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16 

26 

244.8805 

-35.2587 

17 

26 

244.8939 

-36.5104 

18 

26 

244.9278 

-37.5468 

19 

26 

244.9745 

-38.3838 

20 

26 

245.0293 

-39.0368 

1 

27 

241.1034 

-1.4449 

2 

27 

241.1737 

-4.3133 

3 

27 

241.2641 

-7.1231 

4 

27 

241.3932 

-9.8336 

5 

27 

241.5769 

-12.4048 

6 

27 

241.7908 

-14.8563 

7 

27 

241.9846 

-17.2393 

8 

27 

242.1331 

-19.5599 

9 

27 

242.2261 

-21.8111 

10 

27 

242.2633 

-23.9937 

11 

27 

242.2550 

-26.1054 

12 

27 

242.2221 

-28.1278 

13 

27 

242.1796 

-30.0439 

14 

27 

242.1377 

-31.8316 

15 

27 

242.1069 

-33.4747 

16 

27 

242.0938 

-34.9445 

17 

27 

242.1068 

-36.2040 

18 

27 

242.1395 

-37.2518 

19 

27 

242.1845 

-38.1026 

20 

27 

242.2374 

-38.7712 

1 

28 

238.3611 

-1.4287 

2 

28 

238.4285 

-4.2650 

3 

28 

238.5150 

-7.0441 

4 

28 

238.6387 

-9.7267 

5 

28 

238.8146 

-12.2735 

6 

28 

239.0196 

-14.7022 

7 

28 

239.2053 

-17.0610 

8 

28 

239.3476 

-19.3557 

9 

28 

239.4367 

-21.5803 

10 

28 

239.4724 

-23.7377 

11 

28 

239.4644 

-25.8277 

12 

28 

239.4329 

-27.8325 

13 

28 

239.3922 

-29.7356 

14 

28 

239.3520 

-31.5152 

15 

28 

239.3225 

-33.1550 

16 

28 

239.3100 

-34.6266 

17 

28 

239.3224 

-35.8927 

18 

28 

239.3538 

-36.9514 

19 

28 

239.3969 

-37.8161 

20 

28 

239.4476 

-38.5010 

1 

29 

235.6268 

-1.4123 

2 

29 

235.6907 

-4.2165 

3 

29 

235.7730 

-6.9650 

4 

29 

235.8906 

-9.6201 

5 

29 

236.0579 

-12.1433 

6 

29 

236.2528 

-14.5505 

7 

29 

236.4294 

-16.8863 

60.2979 

47.0509 

8.5118 

57.6065 

50.2746 

8.4711 

54.8391 

53.4969 

8.3803 

52.0183 

56.7224 

8.3467 

49.1509 

59.9663 

8.3324 

85.8093 

1.4451 

8.7133 

85.2552 

4.3408 

8.7055 

84.4769 

7.2324 

8.6905 

83.4212 

10.1078 

8.6726 

82.0695 

12.9595 

8.6681 

80.5022 

15.8163 

8.6826 

78.8234 

18.7144 

8.6761 

77.0546 

21.6573 

8.6625 

75.1960 

24.6410 

8.6277 

73.2513 

27.6672 

8.6164 

71.2173 

30.7373 

8.6020 

69.0833 

33.8446 

8.5879 

66.8411 

36.9856 

8.5736 

64.4903 

40.1511 

8.5312 

62.0321 

43.3382 

8.5147 

59.4750 

46.5325 

8.4227 

56.8220 

49.7239 

8.3795 

54.0931 

52.9132 

8.2875 

51.3104 

56.1049 

8.2521 

48.4805 

59.3143 

8.2360 

84.6419 

1.4290 

8.6241 

84.0878 

4.2922 

8.6149 

83.3174 

7.1521 

8.5978 

82.2775 

9.9961 

8.5779 

80.9475 

12.8169 

8.5716 

79.4029 

15.6422 

8.5859 

77.7440 

18.5066 

8.5785 

75.9929 

21.4139 

8.5646 

74.1520 

24.3608 

8.5278 

72.2278 

27.3497 

8.5172 

70.2183 

30.3831 

8.5029 

68.1123 

33.4546 

8;4918 

65.9015 

36.5606 

8.4794 

63.5848 

39.6918 

8.4392 

61.1626 

42.8453 

8.4239 

58.6430 

46.0063 

8.3325 

56.0284 

49.1642 

8.2877 

53.3383 

52.3195 

8.1947 

50.5942 

55.4763 

8.1574 

47.8022 

58.6502 

8.1396 

83.4508 

1.4126 

8.5341 

82.8982 

4.2435 

8.5235 

82.1372 

7.0713 

8.5048 

81.1150 

9.8842 

8.4825 

79.8093 

12.6746 

8.4736 

78.2906 

15.4689 

8.4877 

76.6542 

18.3000 

8.4800 

D-12 


8 

29 

236.5648 

-19.1559 

9 

29 

236.6495 

-21.3543 

10 

29 

236.6834 

-23.4866 

11 

29 

236.6758 

-25.5542 

12 

29 

236.6459 

-27.5401 

13 

29 

236.6071 

-29.4289 

14 

29 

236.5689 

-31.1985 

15 

29 

236.5408 

-32.8334 

16 

29 

236.5290 

-34.3049 

17 

29 

236.5407 

-35.5765 

18 

29 

236.5705 

-36.6456 

19 

29 

236.6116 

-37.5243 

20 

29 

236.6599 

-38.2261 

1 

30 

232.9003 

-1.3958 

2 

30 

232.9606 

-4.1677 

3 

30 

233.0382 

-6.8857 

4 

30 

233.1490 

-9.5138 

5 

30 

233.3067 

-12.0143 

6 

30 

233.4904 

-14.4011 

7 

30 

233.6568 

-16.7152 

8 

30 

233.7845 

-18.9606 

9 

30 

233.8644 

-21.1332 

10 

30 

233.8964 

-23.2403 

11 

30 

233.8893 

-25.2848 

12 

30 

233.8611 

-27.2509 

13 

30 

233.8245 

-29.1238 

14 

30 

233.7885 

-30.8818 

15 

30 

233.7619 

-32.5096 

16 

30 

233.7508 

-33.9795 

17 

30 

233.7619 

-35.2553 

18 

30 

233.7900 

-36.3344 

19 

30 

233.8286 

-37.2272 

20 

30 

233.8742 

-37.9466 

1 

31 

230.1818 

-1.3793 

2 

31 

230.2381 

-4.1186 

3 

31 

230.3104 

-6.8063 

4 

31 

230.4138 

-9.4076 

5 

31 

230.5609 

-11.8864 

6 

31 

230.7324 

-14.2541 

7 

31 

230.8877 

-16.5475 

8 

31 

231.0069 

-18.7696 

9 

31 

231.0815 

-20.9170 

10 

31 

231.1113 

-22.9988 

11 

31 

231.1047 

-25.0196 

12 

31 

231.0783 

-26.9646 

13 

31 

231.0442 

-28.8202 

14 

31 

231.0106 

-30.5648 

15 

31 

230.9858 

-32.1838 

16 

31 

230.9754 

-33.6502 

17 

31 

230.9857 

-34.9294 

18 

31 

231.0120 

-36.0177 

19 

31 

231.0481 

-36.9249 

74.9226 

21.1718 

8.4657 

73.1004 

24.0816 

8.4274 

71.1971 

27.0329 

8.4172 

69.2119 

30.0289 

8.4030 

67.1336 

33.0633 

8.3943 

64.9537 

36.1331 

8.3832 

62.6705 

39.2286 

8.3450 

60.2840 

42.3470 

8.3311 

57.8017 

45.4732 

8.2405 

55.2257 

48.5964 

8.1944 

52.5748 

51.7165 

8.1011 

49.8696 

54.8377 

8.0623 

47.1161 

57.9750 

8.0428 

82.2361 

1.3960 

8.4432 

81.6863 

4.1945 

8.4308 

80.9363 

6.9903 

8.4110 

79.9338 

9.7722 

8.3862 

78.6552 

12.5326 

8.3742 

77.1651 

15.2965 

8.3876 

75.5539 

18.0948 

8.3795 

73.8437 

20.9312 

8.3654 

72.0414 

23.8037 

8.3260 

70.1592 

26.7170 

8.3163 

68.1981 

29.6748 

8.3025 

66.1470 

32.6714 

8.2953 

63.9975 

35.7038 

8.2856 

61.7473 

38.7623 

8.2489 

59.3963 

41.8438 

8.2362 

56.9512 

44.9338 

8.1472 

54.4140 

48.0209 

8.1001 

51.8026 

51.1050 

8.0071 

49.1367 

54.1898 

7.9673 

46.4221 

57.2898 

7.9461 

80.9977 

1.3795 

8.3512 

80.4522 

4.1450 

8.3379 

79.7147 

6.9089 

8.3168 

78.7339 

9.6600 

8.2890 

77.4849 

12.3910 

8.2738 

76.0266 

15.1250 

8.2860 

74.4432 

17.8908 

8.2777 

72.7562 

20.6920 

8.2637 

70.9748 

23.5274 

8.2240 

69.1142 

26.4023 

8.2142 

67.1769 

29.3213 

8.2008 

65.1525 

32.2791 

8.1950 

63.0330 

35.2728 

8.1859 

60.8155 

38.2928 

8.1508 

58.4995 

41.3363 

8.1390 

56.0915 

44.3885 

8.0521 

53.5933 

47.4384 

8.0050 

51.0216 

50.4856 

7.9129 

48.3954 

53.5333 

7.8721 

D-13 


20 

31 

231.0906 

-37.6624 

1 

32 

227.4712 

-1.3626 

2 

32 

227.5231 

-4.0693 

3 

32 

227.5898 

-6.7267 

4 

32 

227.6851 

-9.3017 

5 

32 

227.8207 

-11.7596 

6 

32 

227.9788 

-14.1095 

7 

32 

228.1220 

-16.3834 

8 

32 

228.2319 

-18.5830 

9 

32 

228.3007 

-20.7056 

10 

32 

228.3282 

-22.7621 

11 

32 

228.3221 

-24.7585 

12 

32 

228.2978 

-26.6814 

13 

32 

228.2663 

-28.5181 

14 

32 

228.2353 

-30.2476 

15 

32 

228.2125 

-31.8560 

16 

32 

228.2029 

-33.3173 

17 

32 

228.2124 

-34.5985 

18 

32 

228.2366 

-35.6957 

19 

32 

228.2699 

-36.6171 

20 

32 

228.3091 

-37.3736 

1 

33 

224.7686 

-1.3458 

2 

33 

224.8156 

-4.0198 

3 

33 

224.8762 

-6.6469 

4 

33 

224.9628 

-9.1960 

5 

33 

225.0860 

-11.6340 

6 

33 

225.2295 

-13.9672 

7 

33 

225.3597 

-16.2228 

8 

33 

225.4595 

-18.4008 

9 

33 

225.5220 

-20.4991 

10 

33 

225.5470 

-22.5302 

11 

33 

225.5414 

-24.5016 

12 

33 

225.5194 

-26.4014 

13 

33 

225.4908 

-28.2176 

14 

33 

225.4626 

-29.9303 

15 

33 

225.4419 

-31.5261 

16 

33 

225.4331 

-32.9807 

17 

33 

225.4418 

-34.2627 

18 

33 

225.4638 

-35.3683 

19 

33 

225.4940 

-36.3042 

20 

33 

225.1)297 

-37.0801 

1 

34 

222.0738 

-1.3289 

2 

34 

222.1158 

-3.9701 

3 

34 

222.1697 

-6.5669 

4 

34 

222.2469 

-9.0906 

5 

34 

222.3567 

-11.5094 

6 

34 

222.4847 

-13.8272 

7 

34 

222.6008 

-16.0657 

8 

34 

222.6898 

-18.2231 

9 

34 

222.7455 

-20.2976 

10 

34 

222.7678 

-22.3032 

11 

34 

222.7628 

-24.2488 

45.7203 

56.5951 

7.8492 

79.7357 

1.3629 

8.2585 

79.1959 

4.0953 

8.2443 

78.4724 

6.8272 

8.2218 

77.5152 

9.5477 

8.1915 

76.2985 

12.2496 

8.1724 

74.8751 

14.9544 

8.1829 

73.3220 

17.6883 

8.1745 

71.6602 

20.4544 

8.1607 

69.9007 

23.2526 

8.1213 

68.0619 

26.0889 

8.1116 

66.1482 

28.9686 

8.0983 

64.1503 

31.8865 

8.0927 

62.0602 

34.8404 

8.0843 

59.8749 

37.8206 

8.0504 

57.5937 

40.8245 

8.0394 

55.2227 

43.8376 

7.9555 

52.7634 

46.8493 

7.9087 

50.2318 

49.8586 

7.8183 

47.6458 

52.8683 

7.7769 

45.0106 

55.8914 

7.7530 

78.4501 

1.3459 

8.1647 

77.9174 

4.0454 

8.1498 

77.2093 

6.7452 

8.1265 

76.2777 

9.4352 

8.0939 

75.0961 

12.1086 

8.0707 

73.7104 

14.7847 

8.0790 

72.1904 

17.4870 

8.0703 

70.5555 

20.2185 

8.0567 

68.8191 

22.9794 

8.0183 

67.0024 

25.7769 

8.0081 

65.1122 

28.6166 

7.9949 

63.1401 

31.4941 

7.9887 

61.0790 

34.4069 

7.9809 

58.9254 

37.3461 

7.9482 

56.6788 

40.3087 

7.9375 

54.3447 

43.2814 

7.8571 

51.9245 

46.2536 

7.8112 

49.4333 

49.2242 

7.7237 

46.8878 

52.1952 

7.6817 

44.2930 

55.1789 

7.6568 

77.1408 

1.3291 

8.0699 

76.6166 

3.9953 

8.0548 

75.9255 

6.6629 

8.0306 

75.0215 

9.3226 

7.9958 

73.8776 

11.9678 

7.9686 

72.5327 

14.6158 

7.9741 

71.0483 

17.2870 

7.9649 

69.4423 

19.9841 

7.9512 

67.7299 

22.7079 

7.9143 

65.9358 

25.4662 

7.9030 

64.0688 

28.2653 

7.8903 
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12 

34 

222.7431 

-26.1243 

62.1222 

31.1014 

7.8833 

13 

34 

222.7176 

-27.9186 

60.0895 

33.9722 

7.8756 

14 

34 

222.6925 

-29.6127 

57.9672 

36.8689 

7.8436 

15 

34 

222.6740 

-31.1942 

55.7548 

39.7890 

7.8329 

16 

34 

222.6662 

-32.6402 

53.4574 

42.7197 

7.7569 

17 

34 

222.6740 

-33.9220 

51.0765 

45.6512 

7.7127 

18 

34 

222.6936 

-35.0354 

48.6261 

48.5820 

7.6286 

19 

34 

222.7206 

-35.9859 

46.1215 

51.5137 

7.5869 

20 

34 

222.7523 

-36.7820 

43.5676 

54.4572 

7.5612 

1 

35 

219.3870 

-1.3119 

75.8079 

1.3122 

7.9743 

2 

35 

219.4235 

-3.9201 

75.2936 

3.9448 

7.9587 

3 

35 

219.4704 

-6.4868 

74.6211 

6.5803 

7.9343 

4 

35 

219.5375 

-8.9854 

73.7465 

.9.2098 

7.8981 

5 

35 

219.6330 

-11.3861 

72.6430 

11.8274 

7.8668 

6 

35 

219.7443 

-13.6897 

71.3420 

14.4478 

7.8691 

7 

35 

219.8452 

-15.9122 

69.8957 

17.0882 

7.8585 

8 

35 

219.9227 

-18.0497 

68.3205 

19.7509 

7.8448 

9 

35 

219.9712 

-20.1008 

66.6332 

22.4375 

7.8100 

10 

35 

219.9906 

-22.0809 

64.8619 

25.1566 

7.7973 

11 

35 

219.9862 

-24.0002 

63.0180 

27.9145 

7.7849 

12 

35 

219.9691 

-25.8503 

61.0964 

30.7085 

7.7763 

13 

35 

219.9469 

-27.6212 

59.0916 

33.5360 

7.7682 

14 

35 

219.9250 

-29.2950 

57.0002 

36.3890 

7.7365 

15 

35 

219.9089 

-30.8601 

54.8217 

39.2649 

7.7257 

16 

35 

219.9021 

-32.2960 

52.5610 

42.1522 

7.6549 

17 

35 

219.9089 

-33.5764 

50.2195 

45.0418 

7.6129 

18 

35 

219.9260 

-34.6971 

47.8101 

47.9318 

7.5334 

19 

35 

219.9494 

-35.6623 

45.3469 

50.8233 

7.4924 

20 

35 

219.9771 

-36.4792 

42.8344 

53.7257 

7.4664 

1 

36 

216.7081 

-1.2949 

74.4513 

1.2950 

7.8778 

2 

36 

216.7387 

-3.8700 

73.9483 

3.8941 

7.8624 

3 

36 

216.7781 

-6.4066 

73.2959 

6.4974 

7.8379 

4 

36 

216.8345 

-8.8804 

72.4528 

9.0968 

7.8006 

5 

36 

216.9147 

-11.2638 

71.3923 

11.6871 

7.7663 

6 

36 

217.0082 

-13.5545 

70.1382 

14.2803 

7.7639 

7 

36 

217.0931 

-15.7622 

68.7327 

16.8903 

7.7517 

8 

36 

217.1581 

-17.8808 

67.1901 

19.5188 

7.7374 

9 

36 

217.1989 

-19.9090 

65.5290 

22.1683 

7.7049 

10 

36 

217.2152 

-21.8634 

63.7809 

24.8477 

7.6907 

11 

36 

217.2116 

-23.7558 

61.9599 

27.5640 

7.6788 

12 

36 

217.1972 

-25.5793 

60.0628 

30.3148 

7.6677 

13 

36 

217.1785 

-27.3253 

58.0855 

33.0980 

7.6591 

14 

36 

217.1602 

-28.9772 

56.0244 

35.9061 

7.6274 

15 

36 

217.1466 

-30.5241 

53.8796 

38.7362 

7.6162 

16 

36 

217.1409 

-31.9481 

51.6553 

41.5787 

7.5509 

17 

36 

217.1466 

-33.2259 

49.3534 

44.4249 

7.5116 

18 

36 

217.1609 

-34.3534 

46.9854 

47.2731 

7.4376 

19 

36 

217.1807 

-35.3334 

44.5639 

50.1233 

7.3983 

20 

36 

217.2039 

-36.1718 

42.0933 

52.9836 

7.3725 

1 

37 

214.0371 

-1.2777 

73.0711 

1.2778 

7.7798 

2 

37 

214.0615 

-3.8195 

72.5808 

3.8431 

7.7651 

3 

37 

214.0930 

-6.3261 

71.9499 

6.4140 

7.7410 
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4 

37 

214.1380 

-8.7757 

5 

37 

214.2020 

-11.1426 

6 

37 

214.2766 

-13.4216 

7 

37 

214.3443 

-15.6158 

8 

37 

214.3962 

-17.7162 

9 

37 

214.4288 

-19.7220 

10 

37 

214.4418 

-21.6508 

11 

37 

214.4389 

-23.5155 

12 

37 

214.4274 

-25.3114 

13 

37 

214.4125 

-27.0310 

14 

37 

214.3979 

-28.6591 

15 

37 

214.3871 

-30.1861 

16 

37 

214.3825 

-31.5964 

17 

37 

214.3871 

-32.8705 

18 

37 

214.3985 

-34.0044 

19 

37 

214.4142 

-34.9992 

20 

37 

214.4328 

-35.8597 

1 

38 

211.3740 

-1.2604 

2 

38 

211.3919 

-3.7689 

3 

38 

211.4149 

-6.2456 

4 

38 

211.4479 

-8.6712 

5 

38 

211.4948 

-11.0226 

6 

38 

211.5494 

-13.2911 

7 

38 

211.5990 

-15.4728 

8 

38 

211.6370 

-17.5560 

9 

38 

211.6608 

-19.5399 

10 

38 

211.6704 

-21.4429 

11 

38 

211.6683 

-23.2794 

12 

38 

211.6598 

-25.0466 

13 

38 

211.6489 

-26.7383 

14 

38 

211.6382 

-28.3409 

15 

38 

211.6303 

-29.8459 

16 

38 

211.6269 

-31.2409 

17 

38 

211.6302 

-32.5103 

18 

38 

211.6386 

-33.6499 

19 

38 

211.6502 

-34.6597 

20 

38 

211.6638 

-35.5429 

1 

39 

208.7188 

-1.2429 

2 

39 

208.7298 

-3.7180 

3 

39 

208.7440 

-6.1649 

4 

39 

208.7642 

-8.5669 

5 

39 

208.7930 

-10.9038 

6 

39 

208.8266 

-13.1630 

7 

39 

208.8570 

-15.3334 

8 

39 

208.8804 

-17.4002 

9 

39 

208.8950 

-19.3627 

10 

39 

208.9009 

-21.2399 

11 

39 

208.8996 

-23.0474 

12 

39 

208.8944 

-24.7848 

13 

39 

208.8877 

-26.4471 

14 

39 

208.8811 

-28.0224 

15 

39 

208.8763 

-29.5037 

71.1403 

8.9836 

7.7038 

70.1256 

11.5468 

7.6667 

68.9213 

14.1131 

7.6593 

67.5592 

16.6931 

7.6442 

66.0511 

19.2874 

7.6288 

64.4173 

21.8996 

7.5991 

62.6926 

24.5392 

7.5832 

60.8943 

27.2132 

7.5715 

59.0214 

29.9201 

7.5574 

57.0710 

32.6577 

7.5477 

55.0399 

35.4194 

7.5159 

52.9283 

38.2025 

7.5038 

50.7405 

40.9983 

7.4449 

48.4782 

43.7999 

7.4097 

46.1519 

46.6049 

7.3419 

43.7726 

49.4126 

7.3046 

41.3443 

52.2297 

7.2797 

71.6672 

1.2607 

7.6812 

71.1911 

3.7919 

7.6678 

70.5833 

6.3305 

7.6448 

69.8091 

8.8701 

7.6081 

68.8428 

11.4065 

7.5693 

67.6913 

13.9463 

7.5559 

66.3753 

16.4960 

7.5364 

64.9035 

19.0562 

7.5196 

63.2980 

21.6312 

7.4924 

61.5972 

24.2304 

7.4744 

59.8214 

26.8616 

7.4635 

57.9721 

29.5237 

7.4458 

56.0483 

32.2144 

7.4345 

54.0466 

34.9283 

7.4025 

51.9680 

37.6625 

7.3890 

49.8165 

40.4101 

7.3370 

47.5940 

43.1656 

7.3061 

45.3096 

45.9260 

7.2461 

42.9729 

48.6899 

7.2114 

40.5875 

51.4622 

7.1881 

70.2397 

1.2432 

7.5815 

69.7792 

3.7403 

7.5696 

69.1960 

6.2465 

7.5481 

68.4592 

8.7563 

7.5135 

67.5439 

11.2661 

7.4747 

66.4483 

13.7794 

7.4547 

65.1809 

16.2990 

7.4290 

63.7474 

18.8248 

7.4099 

62.1712 

21.3622 

7.3857 

60.4945 

23.9208 

7.3649 

58.7410 

26.5083 

7.3545 

56.9150 

29.1246 

7.3320 

55.0172 

31.7673 

7.3192 

53.0445 

34.4319 

7.2865 

50.9987 

37.1156 

7.2714 
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16 

39 

208.8742 

-30.8818 

48.8833 

39.8134 

17 

39 

208.8762 

-32.1451 

46.7008 

42.5208 

18 

39 

208.8813 

-33.2900 

44.4586 

45.2351 

19 

39 

208.8884 

-34.3148 

42.1649 

47.9535 

20 

39 

208.8968 

-35.2215 

39.8228 

50.6795 

1 

40 

206.0716 

-1.2254 

68.7885 

1.2257 

2 

40 

206.0753 

-3.6668 

68.3450 

3.6884 

3 

40 

206.0802 

-6.0840 

67.7880 

6.1621 

4 

40 

206.0870 

-8.4629 

67.0904 

8.6421 

5 

40 

206.0968 

-10.7860 

66.2289 

11.1252 

6 

40 

206.1081 

-13.0372 

65.1922 

13.6120 

7 

40 

206.1185 

-15.1975 

63.9760 

16.1013 

8 

40 

206.1263 

-17.2489 

62.5826 

18.5925 

9 

40 

206.1313 

-19.1904 

61.0369 

21.0919 

10 

40 

206.1333 

-21.0417 

59.3847 

23.6093 

11 

40 

206.1329 

-22.8196 

57.6532 

26.1525 

12 

40 

206.1311 

-24.5260 

55.8501 

28.7218 

13 

40 

206.1288 

-26.1574 

53.9777 

31.3153 

14 

40 

206.1266 

-27.7038 

52.0336 

33.9291 

15 

40 

206.1250 

-29.1595 

50.0202 

36.5605 

16 

40 

206.1243 

-30.5189 

47.9409 

39.2066 

17 

40 

206.1250 

-31.7751 

45.7985 

41.8642 

18 

40 

206.1267 

-32.9247 

43.5988 

44.5304 

19 

40 

206.1291 

-33.9647 

41.3486 

47.2015 

20 

40 

206.1319 

-34.8955 

39.0503 

49.8791 

7.2268 

7.2009 

7.1501 

7.1189 

7.0982 

7.4807 

7.4710 

7.4525 

7.4212 

7.3838 

7.3557 

7.3220 

7.2992 

7.2781 

7.2542 

7.2447 

7.2167 

7.2023 

7.1679 

7.1509 

7.1149 

7.0949 

7.0541 

7.0272 

7.0096 
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APPENDIX  E 

COMPUTER  CODE  INTERP 


PROGRAM  INTERP 
$DEBUG 


This  program  was  written  to  support  the  B-1B  Transparency  Test, 
JON  24010505,  May  1990. 


For  a  given  Mach  number,  altitude,  and  angle  cf  attack,  this  program 
interpolates  Cp  data  from  Mach  number  and  alpha  tables  to  produce 
pressure  in  psi  at  panel  centroids. 

This  is  version  2.  It  creates  plot  files  for  TECPLOT  internally. 

REAL  MACH, MACHO, MAX, MIN, MSTAR 
CHARACTER* 4 6  TITLE 

CHARACTER* 1 0  FILENAME ( 6 ) , TAPE , FILE2 , FILE3 , FILE4 , FILE5 
DIMENSION  CP(20,40) ,DP(20,40) ,DT(20,40) , TEMP (40) 

DIMENSION  ROATM(lOl) , PATM(lOl) , SSATM( 101) , TATM(lOl) ,A(6) 

DIMENSION  MACH ( 6 ) , N ALPHA ( 6 ) , ALPHA (6,5), COUNT ( 6 ) 

NTABLE-6 

GAMMA-1.4 

EXP- (GAMMA-1 . ) /GAMMA 
DATA  CP/800*0.0/ 

DATA  MACH/. 30, .56, .70, .75, .  ',.91/ 

DATA  NALPHA/ 3, 2, 4, 3, 3, 3/ 

DATA  ALPHA/0.,0.,0.,0.,0.,0., 

*  4. ,2. , 1. , 1. , 1. , 1. , 

*  S.,0.,2.,2.,2.,2., 

*  0. ,0. ,3. ,0. ,0. ,0. , 

*  0. ,  0 . , 0 « ,  0 . ,  0 . ,  0 .  / 

DATA  COUNT/ 1. ,2. ,3. ,4. ,5. ,6,/ 

Atmospheric  data 

DATA  ROATM  /  .23770,-23080, 

*.22410,-21750, .21110, .20480, . 19870, . 19270, . 18680, . 18110, 

* . 17  550 , . 17010 , . 16480 ,. 15960 ,. 15450 ,. 14960 , . 14470 , . 14000 , 

* . 13550, . 13100, . 12660, . 12240, . 11830, .11420, . 11030, . 10650, 

*.10280, .09920, .09570, .09230, .08890, .08570, .08260, .07950,. 

*.07650, .07370, .07090, .06760,-06440, .06140, .05850, .05580, 

* . 05310 , . 05070 , . 04830 , . 04600, . 04390 , . 04180, . 03980, . 03800, 

*.03620, .03448, .03286, .03132, .02985, .02845, .02711, .02584, 

*.02463, .02347, .02237, .02132, .02032, .01937, .01846, .01759, 

*.01676, .01595, .01518, .01445, .01376, .01310, .01247, .01187, 

*.01131, .01077, .01025, .00976, .00930, .00886, .00844, .00804, 

*.00767, .00730, .00696,-00663, .00632, .00602, .00574, .00547, 

* . 00522 , . 00497 ,. 00474 ,. 00452 ,. 0043 1 ,. 0041 1 , . 00392 , . 00374 , 

*.00357, .00340, .00324/ 

DATA  PATH  /  2116.2,2040.8, 


*1967.7,1896.6, 

1827.7, 

1760.8, 

1695.9, 

1632.9, 

1571.9, 

1512.7 

*1455.3,1399.7, 

1345.9, 

1293.7, 

1243.2, 

1194.3, 

1146.9, 

1101.1 

*1056.8,1013.9, 

972.5, 

932.4, 

893.7, 

856.3, 

820.2, 

785.3 

*  751.6,  719.1, 

687.8, 

657.6, 

628.4, 

600.3, 

573.3, 

547.2 

•  522.1,  498.0, 

474.7, 

452.4, 

431.2, 

411.0, 

391.7, 

373.3 

*  355.8,  339.1, 

323.2, 

308.0, 

293.5, 

279.8, 

266.6, 

254.1 

•  242.2,  230.8, 

220.0, 

209.7, 

199.9, 

190.5, 

181.5, 

173.0 

*  164.9,  157.2, 

149.8, 

142.8, 

136.1, 

129.7, 

123.6, 

117.8 
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*  112. 

3,  107.0,  102 

.0,  97. 

2,  92 

.7,  88.4,  84.3,  80. 

3, 

* 

76.6, 

73.1, 

69.7, 

66.5, 

63.4,  60.5,  57.7, 

55.0 

* 

52.5, 

50.1, 

47.8, 

45.6, 

43.5,  41.5,  39.6, 

37.8 

* 

* 

DATA 

36.1,  34.5, 

25.0,  23.8, 

SSATM  / 

32.9, 

22.8/ 

31.4, 

30.0,  28.7,  27.4, 

661.03,659.20, 

26.1 

*656.92,654.63,652.32,650.00,647.69,645.36,643.02,640.69, 

*638.34,635.97,633.60,631.22,628.83,626.44,624.03,621.62, 

*619.19,616.76,614.31,611.86,609.40,606.93,604.44,601.95, 

*599.44,596.93,594.41,591.86,589.32,586.77,584.20,581.62, 

*579.03,576.43,573.80,573.57,573.57,573.57,573.57,573.57, 

*573.57,573.57,573.57,573.57,573.57,573.57,573.57,573.57, 

*573.57,573.16,573.16,573.16,573.16,573.16,573.16,573.16, 

*573.16,573.16,573.16,573.16,573.16,573.16,573.16,573.16, 

*573.32,573.72,574.12,574.53,574.92,575.33,575.73,576.13, 

*576.53,576.94,577.37,577.73,578.14,578.53,578.93,579.33, 

*579.73,580.13,580.53,580.93,581.33,581.72,582.12,582.52, 

*582.91,583.31,583.71,584.10,584.50,584.89,585.29,585.68, 

*586.08,586.47,586.86/ 

DATA  TATM  /  518.67,515.10, 

*511.54,507.97,504.40,500.84,497.27,493.71,490.14,486.57, 
*483.01,479.44,475.87,472.31,468.74,465.18,461.61,458.04, 
*454.48,450.91,447.35,443.78,440.21,436.65,433.08,429.52, 
*425.95,422.38,418.82,415.25,411.68,408.12,404.55,400.99, 
*397.41,393.85,390.29,389.97,389.97,389.97,389.97,389.97, 
*389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 
*389.97,389.97,389.97,389.97,389.97,389.97,389.97,389.97, 
*389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 . 97 , 389 .97 , 389 . 97 , 389 . 97 , 
*390.18,390.73,391.28,391.83,392.37,392.92,393.47,394.02, 
*394.57,395.12,395.67,396.22,396.76,397.31,397.86,398.41, 
*398.96,399.51,400.06,400.60,401.15,401.70,402.30,402.80, 
*403.35,403.90,404.44,404.99,405.54,406.90,406.64,407.19, 
*407.74,408.28,408.83/ 

Define  input  files  containing  areo  data 

FILENAME ( 1 ) = ' CANOPY . M30 ' 

FILENAMB ( 2 ) - ' CANOPY .  M56 ' 

FILENAME ( 3 ) * ' CANOPY . M70 ' 

FILENAME ( 4 )  * ' CANOPY . M75 ' 

FILENAME ( 5 ) « ' CANOPY . M85 ' 

FILENAME  ( 6 ) « '  CANOPY .  M9 1 ' 

Read  the  flight  condition 

WRITE (*,*)*  INPUT  MACH  NUMBER  (F10)' 

READ( *, ' (F10.3) ' ) MACHO 

WRITE (*,*)'  INPUT  ANGLE  OF  ATTACK  OR  GROSS  WEIGHT/1000' 

READ ( * , ' (F10.3) ' )ALPHA0 

IF ( ALPHAO . LE .30.0)  GOTO  2 

WRITE (*,*)'  INPUT  WING  SWEEP  ANGLE  (F10) ' 

READ ( * , ' (F10.3) ' ) SWEEPLE 
2  CONTINUE 

WRITE (*,*)'  INPUT  ALTITUDE  DIVIDED  BY  1000' 

READ( *, ' (F10.3) ' ) ALTO 

Obtain  properties  of  the  atmosphere  at  altitude 

NALT»IFIX(ALT0)+1 
RHO-ROATM { NALT ) / 100 . 

PO-PATM(NALT) 

PINT-PO 
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IF(PINT.LT. 1571.877)  PINT-1571 . 877 
SS-SSATM(NALT) *1.689 
TO-TATM { NALT ) 

Q=0. 5*RH0* (MACHO*SS )  **2 
WRITE{*, ' (4H  Q  F10. 5 , 4H  PSF) ' )Q 

Calculate  angle  of  attack  if  not  supplied 

IF(ALPHAO.LE. 30.0)  GOTO  3 
WEIGHT-1000 . 0* ALPHA 0 

AREA-1907 . 0+ ( 22 59 . 0-1907 . 0 )*( SWEEPLE-15 . 0 )/( 65 . -15 . ) 

SPAN— (254. 0+1386. 0*COS( (SWEEPLE-6.2)/57.2958) )/12. 

AR=SP AN* SPAN /AREA 
SWEEP— SWEEPLE- 5 . 5 

TPR=0.35+( .526-0. 35)* (SWEEPLE-15. 0)/( 65.-15. ) 

CALL  CLALPHA( MACHO, AR, SWEEP , TPR , CLA , MSTAR ) 

CL-WEIGHT/Q/AREA 

Alpha  Body  -  (CL/CLA)  -  (Alpha  zero  lift  of  wing)  -  (Alpha  wing/body) 
Alpha  zero  lift  of  wing  may  be  a  weak  function  of  Mach  number 
Assume : 

ALP ZERO  -  0.5 

Alpha  wing/body  may  be  a  function  of  wing  sweep 
Assume : 

ALPWB  -  1.5 

ALPHA0— CL/CLA*57.2958  -ALP ZERO  -  ALPWB 

WRITE ( * »  * ) '  MACH, WEIGHT, SWEEP, ALPHA, CL, CLA, AREA, SPAN, AR, TPR' 

WRITE ( * , 6000 ) MACHO , WEIGHT , SWEEPLE , ALPHA0 , CL , CLA , AREA , SPAN , AR , TPR 
6000  FORMAT (F6.3,F8.0,F5.1,F6.2,F6.3,F6.3,F7.0,F6.1,F6.3,F5.3) 

3  CONTINUE 

Find  the  closest  Mach  number  and  name  of  the  datafile. 

MIN-100 

DO  10  1*1, NT ABLE 
TEST-ABS ( MACHO -MACH ( I ) ) 

IF ( TEST. CT. MIN)  GOTO  10 
IMACH-I 
MIN-TEST 
10  CONTINUE 

Select  pair  of  alpha  tables  for  interpolation/extrapolation 

N-NALPHA ( I MACH ) 

DO  5  1*1, N 

5  A ( I ) -ALPHA ( I MACH , I ) 

CALL  INTER (A, COUNT, N , ALPHA0 , FACTR ) 

Calculate  interpolation/extrapolation  multipliers  for  pair  of  tables 
NTABL1-IFIX ( FACTR) 
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NTABL2 *NTABL 1+1 
FACTR2-FACTR- FLOAT ( NTABL1 ) 

IF(NTABL2.GT.N)  NX-NTABL2-N 
IF ( NTABLl. LT.l)  NX=NTABL1-1 

IF(NX.NE.O)WRITE(*,*) '  CAUTION  -  EXTRAPOLATING  DATA* 

FACTR2  * FACTR2 +NX 
FACTR1=1 . -FACTR2 
NTABL1*NTABL1“NX 
NTABL2  =NTABL2 -NX 
C 

I=IMACH  „  .  „  . 

WRITE ( * , 1002 ) ALPHA ( I , NTABLl ) , FACTR1 , ALPHA  ( I , NTABL2 ) , FACTR2 
1002  FORMAT ( F3. 0 , 16H  deg  data  tiroes  , F8.4.5H  plus, 

*  F3.0, 16H  deg  data  tiroes  ,F8.4) 

C 

C  Read  the  aero  file  and  scale  the  data 
C 

TAPE-FILENAME ( IMACH ) 

FILE2* ' CP .OUT ' 

FILE3='DP.OUT' 

FILE4*'DT.OUT' 

FILE5= ' MAXMIN . OUT ' 

OPEN ( 1 , FILE=TAPE , STATUS* ' OLD ' ) 

OPEN ( 2 , FILE=FILE2 , STATUS* ' NEW ' ) 

OPEN ( 3 , FILE*FILE3 , STATUS* ' NEW ' ) 

OPEN ( 4, FILE-FILE4, STATUS- 'NEW' ) 

OPEN ( 5 , FILE-FILE5 , STATUS* ' NEW ' ) 

DO  20  K=1 , 6 

IF (K.GT.NTABL2 )  GOTO  30 
FACTR-0.0 

IF ( K . EQ. NTABLl ) FACTR-FACTR1 
IF ( K . EQ. NTABL2 ) FACTR*FACTR2 
READ (1,2000) TITLE 
WRITE ( * , 2000 ) TITLE 
2000  FORMAT(A46) 

DO  20  1*1,20 

READ ( 1 , 2000 ) TITLE 

READ (1,3000) ( TEMP (J),J*1,4Q) 

3000  FORMAT (F9. 6, 7F10. 6) 

3001  FORMAT (F9.5,7F10.5) 

3002  FORMAT(F9.3,7F10.3) 

DO  20  J*l, 40 

20  CP(I, J)*CP(I, J)+TEMP(J)*FACTR 
30  CONTINUE 
C 

C  Calculate  pressure  and  temperature  differential 
C 

C  DP  is  the  inside  pressure  minus  the  outside  pressure  in  PSI 
C 

C  DT  is  the  outside  temperature  in  degrees  F 
C 

DO  35  1-1,20 
DO  35  J-1,40 
CPIJ-CP(I,J) 

DT(I,J)*T0*( (P0+Q*CPIJ)/P0)**EXP  -  459.6 
35  DP ( I , J ) * ( PINT-PO-CPI J*Q ) / 144 . 

C 

C  Output  CP,  DP,  and  DT  on  separate  files 
C 

WRITE ( 2 , 4000 ) MACHO , ALPHAO 
WRITE (3,4001) MACHO , ALPHAO 
WRITE ( 4 , 4002 ) MACHO , ALPHAO 


4000  FORMAT ( 'B1  CANOPY  CP  FOR  MACH  ='F6.3'  AND  ALPHA  *'F5.2) 

4001  FORMAT ( '81  CANOPY  DP  FOR  MACH  ='F6.3'  AND  ALPHA  -'F5.2) 

4002  FORMAT ( 'B1  CANOPY  DT  FOR  MACH  ='F6.3'  AND  ALPHA  *'F5.2) 

C 

DO  40  1*1,20 
WRITE ( 2 , 5000 ) I 
WRITE (3, 5001) I 
WRITE (4, 5002) I 

5000  FORMAT ( ' CP ( I , J )  FOR  I  =  '13'  AND  J  =  1  THROUGH  40') 

5001  FORMAT { ' DP ( I , J )  FOR  I  =  '13'  AND  J  =  1  THROUGH  40') 

5002  FORMAT ( ' DT ( I , «J )  FOR  I  *  '13'  AND  J  *  1  THROUGH  40') 

WRITE (2, 3000) <CP( I, J) , J=l,40) 

WRITE (3, 3001) (DP(I, J) , J=l,40) 

40  WRITE (4, 3002) (DT(I,J),J*1,40) 

C 

C  Locate  max  and  min  of  CP,  DP  and  DT 
C 

WRITE ( 5 , 4000 ) MACHO , ALPHA0 

CALL  MAXMIN(CP,MAX, IMAX, UMAX, MIN, IMIN,JMIN) 

WRITE (*, 5500 ) IMAX , JMAX , MAX , IMIN , JMIN , MIN 
WRITE ( 5 , 5 500 ) IMAX , JMAX , MAX , IMIN , JMIN , MIN 

5500  FORMAT ( '  MAX  IS  CP( ' 12 ' , ' 12 ' ) = ' F10 . 6 

1  ',  MIN  IS  CP ( ' 12 ' , ' 12 ' ) = ' F10. 6// ) 

C 

WRITE ( 5 , 4001 ) MACHO , ALPHAO 

CALL  MAXMIN(DP,MAX, IMAX, JMAX, MIN, IMIN, JMIN) 

WRITE ( * , 5501 ) IMAX, JMAX , MAX, IMIN , JMIN , MIN 
WRITE ( 5 , 5501 ) IMAX, JMAX , MAX , IMIN, JMIN , MIN 

5501  FORMAT ( '  MAX  IS  DP( ' 12 ' , ' 12 ' ) * ' F10 . 6 

1  ',  MIN  IS  DP ( ' 12 ' , ' 12 ' )= ' F10. 6// ) 

C 

WRITE ( 5 , 4002 ) MACHO , ALPHAO 

CALL  MAXMIN(DT,MAX, IMAX, JMAX, MIN, IMIN, JMIN) 

WRITE (*, 5502 ) IMAX, JMAX , MAX , IMIN , JMIN , MIN 
WRITE ( 5 , 5 502 ) IMAX , JMAX , MAX , IMIN , JMIN , MIN 

5502  FORMAT ( '  MAX  IS  DT( ' 12 ' , ' 12 ' )- 'F10. 6 

1  ',  MIN  IS  DT('I2','I2')-'F10.6//) 

C 

C  Create  plot  files  for  CP,  DP,  and  DT 
C 

CLOSE (2) 

CLOSE ( 3 ) 

CLOSE (4) 

CALL  PLOT2 ( FILE2 ) 

CALL  PLOT2 ( FILE3 ) 

CALL  PLOT2 ( FILE4 ) 

STOP 

END 

C 

C 

C 

SUBROUTINE  MAXMIN (A, MAX , IMAX , JMAX , MIN , IMIN , JMIN ) 

C 

C  This  subroutine  locates  the  maximum  and  minimum  values  in 
C  the  array  A  together  with  the  index  of  each. 

C 

REAL  MAX, MIN 
DIMENSION  A(20, 40) 

C 

MAX* A (1,1) 

MIN*A(1,1) 

I MAX-1 
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JMAX-1 

IMIN-1 

JMIN-1 

C 

DO  10  1-1,20 
DO  10  J-l, 40 
C 

IF ( A ( I , J ) . LT . MAX )  GOTO  5 
MAX— A ( I , J ) 

I MAX -I 
JMAX-J 
5  CONTINUE 
C 

IF(A(I,J) . GT.MIN)  GOTO  10 
MIN=A( I , J) 

IMIN-I 
JMIN-J 
10  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  INTER ( X,  Y , NPTS , XVAL , YANS ) 

TWO  DIMENSIONAL  TABLE  LOOKUP 
X-INDEPENDENT  TABLE 
Y-DEPENDENT  TABLE 

NPTS-NUMBER  OF  VALUES  IN  EITHER  X  OR  Y  TABLE 
XVAL- INDEPENDENT  ARGUMENT 
YANS-CALCULATED  RESULT 


DIMENSION  X(l) ,Y(1) 

I -NPTS 

IF(XVAL-X<I) )5,30,35 
IF(XVAL-X(1) )7, 15,20 
IF(X(I)-X(1))8,8,9 
I— NPTS— 1 
GO  TO  50 


9  1*1 

GO  TO  50 
15  YANS-Y(l) 

GO  TO  55 

20  IF(XVAL-X(I) ) 2 5, 2 5, 50 
25  I-I-l 

GO  TO  20 
30  YANS-Y(I) 

GO  TO  55 

35  IF(XVAL-X(1) ) 40, 15, 36 

36  IF(X(I)-X(1))9,9,8 

40  IF(XVAL-X(I) ) 50, 45, 45 
45  I-I-l 

GO  TO  40 

50  YANS— (Y(I+1 )-Y(I) )/ (X(I+1)-X(I) ) * (XVAL-X( I ) )+Y ( I ) 

55  RETURN 
END 


SUBROUTINE  PLOT2 (AEROFILE) 

Thia  varaion  ia  adapt ad  to  TECPLOT  to  produca  contour  plota. 
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C  This  program  reads  the  geometry  file  and  an  aero  file  and  writes  a 
C  file  to  produce  contour  plots  of  Cp,  delta  P,  or  delta  T  on  the 
C  surface  of  the  canopy. 

C 

CHARACTER* 10  TITLE1 
CHARACTER*  2  TITLE2 
CHARACTER* 60  TITLE3 
CHARACTER* 10  AEROFILE 
REAL  MACH 

DIMENSION  CP(20,40) ,X(20,40) ,S(20,40) 

C 

WRITE ( * , 500 ) AEROFILE 

500  FORMAT ( '  Generating  plot  file  for  'A6) 

OPEN ( 1 , FILE-AEROFILE , STATUS* 'OLD ' > 

OPEN ( 2 , FILE* ' CANOPY . SRF ' , STATUS* ' OLD ' ) 

C 

C  Read  the  geometry  file 
C 

READ ( 2 , 1000 ) TITLE1 , TITLE2 , TITLE3 
1000  FORMAT (A10,A2,A60) 

WRITE ( * , 1000 ) TITLE 1, TITLE2 , TITLE3 
READ ( 2 , 1000 ) TITLE3 
C 

DO  10  J*l,40 
DO  10  1*1,20 

READ (2, 2000} II, JJ,XX, YY, ZZ,SS 
2000  FORMAT ( 13, 14, 4F10.4 ) 

X(II,JJ)-XX 
10  S(II,JJ)*SS 
C 

C  Read  the  aero  file,  one  Mach  and  alpha  at  a  time 
C 

WRITE ( • , 3500 ) AEROFILE 

3500  FORMAT( '  Ready  to  read  AEROFILE  as',A10) 

20  READ (1, 3000, END-999 ) TITLE2 , MACH, ALPHA 
3000  FORMAT ( 10X, A2 , 11X, F6 . 3 , 12X, F5 . 2 ) 

WRITE ( • , 4000 ) TITLE2 , MACH , ALPHA 

4000  FORMAT ( '  READING  'A2'  DATA:  MACH  *'F6.3';  ALPHA  *'F5.2) 

C 

C  Read  CP,  DP,  or  DT 
C 

DO  30  1*1,20 
READ(1, 6000)11 
6000  FORMAT(15X,I4) 

30  READ(1,7000) (CP(II, J) , J*l,40) 

7000  FORMAT (F9.6,7F10.6) 

C 

C  Write  plot  file  title  and  variables 
C 

WRITE ( * , 7500 ) TITLE2 
75C0  FORMAT ( '  Writing  plot  file  *A2) 

OPEN ( 3 , PXLE*TITLE2 , STATUS* ' NEW ' ) 

WRITE ( 3 , 8000 ) TITLE2 , MACH , ALPHA , TITLE2 
8000  FORMAT (9HTITLE  -  “,A2,9H  FOR  MACH,F4.2,5H  AND  ,F4.2,5H  DEG”/ 

1  18HVARIABLES  *  X,  S,  A2) 

C 

C  Write  plot  zone  title  and  indices 
C 

WRITE (3, 9000) 

9000  FORMAT (29HZONE  T  *  "ZONE-1”,  1*20,  J*18) 
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OO  40  J*3,38,2 
DO  40  1*1,20 

40  WRITE ( 3 , 10000 ) X ( I , J ) , S ( I , J ) , CP ( I , J ) 

10000  FORMAT (F10.4,F10.4,F10.5) 

GOTO  20 

999  CLOSE ( 3 ) 

RETURN 

END 

C 

SUBROUTINE  CLALPHA  (MACH, AAR, XSWEP, TTPR, X, XMSTR) 

C 

C  This  subroutine  calculates  the  lift  curve  slope  using  an  impirical 
C  formula  used  for  preliminary  design. 

C 

C  AAR-ASPECT  RATIO 

C  MACH-MACH  NUMBER 

C  TTPR-TAPER  RATIO=WING  TIP  CHORD  DIVIDED  BY  ROOT  CHORD 
C  X-LIFT  CURVE  SLOPE 

C  XSWEP-ANGLE  OF  SWEEP  BACK  OF  THE  50  PERCENT  CHORD  LINE 

C 

REAL  MACH 

E  -  8. 0/3.0 

F  =  4. 0/3.0 

G  *  2. 0/3.0 

SSWEP  *  XSWEP/57. 29578 

XMSO  *  (10.  +  .91  *  AAR  **  3) / ( 10.  +  AAR  **  3) 

XMSTR  *  XMSO  +  (1.  -  XMSO)  *  (1.  -  COS(SSWEP))  **  2 

CONE  =  (3.14159  *  AAR)/(1.  +  SQRT{1.  +  ABS(1.  -  (COS(SSWEP)) 

1  **  F  *  (MACH/XMSTR)  **  E)  *  (AAR/ (2.  *  COS(SSWEP)))  **  2)) 

IF  ( MACH. GT. XMSTR)  GO  TO  100 
X  *  CONE 
GOTO  1000 

100  Z  -  ( XMSTR*CONE )  +  ( ( AAR**2 ) / ( ( ( 3 . *3. 14159*AAR) /CONE) * 

1  ( (3 . 14159* AAR) / (CONE)-l. ) * (COS (SSWEP) **G) ) ) 

Y  »  ((1.  +  3 . 14159*AAR) / (3.  +  3.14159  *  AAR))  *  (2.-G  * 

1  SQRT(TTPR)  -  (TTPR)  **  2) 

BE  «  (MACH  -  XMSTR)  *  ((1.  +  (XMSTR /MACH)  **  Y) )  **  2 
X  *  1. /((I. /CONE)  *  (XMSTR/MACH)  **  Z  +  (BE/4)) 

1000  RETURN 
END 
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COMPUTER  CODE  SURF 


PROGRAM  SURF 
C 

C  This  program  determines  the  distance  along  each  cross  section 
C  from  the  top  of  the  canopy  at  the  center  line  to  the  centroid 
C  of  the  panel. 

C 

C  it  is  one  step  in  producing  a  contour  plot  of  pressure  and 
C  temperature  on  the  surface  of  the  canopy, 
c 

CHARACTER*  7 2  TITLE 
REAL  M,M12 
C 

DIMENSION  X(20,40),Y(20,40),Z(20,40),A(20,40) 

C 

OPEN ( 1 , FILE- ' D : CANOPY . GEO ' , STATUS- ' OLD ' ) 

OPEN ( 2 , FILE- ' D : CANOPY . SRF ' , STATUS- ' NEW ' ) 

OPEN ( 3 , FILE- ' D : CANOPY . RAD ' , STATUS- ' NEW ' ) 

WRITE(3,*)'  CURVATURE  DATA:  X,YBAR, ZBAR, RADIUS ' 

READ (1,1000) TITLE 
1000  FORMAT (A72) 

WRITE ( 2 , 1000 ) TITLE 
READ (1,1000) TITLE 

WRITE(2, *) '  I,  J,  XBAR,  YBAR,  ZBAR,  S,  AREA' 

C 

DO  10  1-1,20 
DO  10  J-1,40 

READ ( 1 , 2 000 ) I I , J J , XX , Y Y , Z Z , XCOS , YCOS , ZCOS , AREA 
2000  FORMAT( 13, 14, 3F10. 4, 4F10. 6) 

X(II,JJ)-XX 
Y (II, JJ)-YY 
2(11, JJJ-ZZ 
A(II, JJ)-AREA 
10  CONTINUE 
C 

DO  20  J-1,40 
SUM-0. 

THETAO-O . 0 
C 

C  Calculate  radius  of  curvature  from  three  points  on  a  circle 
C 

C  Get  three  points 
C 

Y1-Y(1,J) 

Z1«Z(1,J) 

Y2— Y(20, J) 

Z2— Z(20, J) 

XX-X(10,J) 

Y3-Y(10, J) 

Z3«Z(10,J) 

C 

C  Find  point  on,  and  slope  of  line  containing  center  of  curvature 
C 

YBAR«(Yl+Y2)/2. 

ZBAR- ( Z1+Z2 ) /2 . 

M12-(Z2-Z1)/ (Y2-Y1) 

M—  1./M12 
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B«ZBAR-M*YBAR 
Find  center  of  curvature 

Y0*.5*(Z3*Z3-Z1*Z1+Y3*Y3-Y1*Y1-2.*B*(Z3-Z1) ) / (M* ( Z3-Z1 ) +Y3-Y1 ) 
ZO*M*YO+B 

RO*SQRT( (Z0-Z1)**2+(Y0-Y1)**2) 

WRITE (3, 2001) XX, YO, Z0,R0 
2001  FORMAT (4F10. 4) 

DO  20  1*1,20 

XX*X(I,J) 

YY*Y ( I, J) 

ZZ*Z(I,J) 

ARE A* A ( I, J) 


R*SQRT( (ZZ-ZBAR) **2+YY*YY) 
THETA=ACOS ( ( ZZ-ZBAR) /R) 
DTHETA«THETA-THETA0 
IF { I . EQ. 1 ) R0*R 
RAVG* ( R+R0 ) / 2 . 

SUM*  SUM+RAVG  *  DTHET A 

SS*SUM 

THETAO*THETA 

WRITE (2 ,2003 ) I , J,XX, YY , ZZ, SS, AREA 
2003  FORMAT (13, 14, 5F10.4) 

20  CONTINUE 
STOP 
END 
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APPENDIX  G 

COMPUTER  CODE  CANOPY 


PROGRAM  CANOPY 
C 

C  This  program  reads  the  QUADPAN  dump  file  and  writes  selected  parameters 
C  to  output  files. 

C 

C  The  X,Y,Z  locations  of  the  centroids  of  the  panels,  the  direction  cosines, 
C  and  the  areas  are  read. 

C 

C  The  pressure  coefficients  are  read  for  each  combination  of  Mach  number  and 
C  angle  of  attack  in  the  QUADPAN  run. 

C 

C  This  program  was  adapted  from  the  program  LOADIT  written  by  John  Riechers 
C  for  use  in  determining  the  pressures  on  the  B-1B  Canopy  for  a  variety  of 
C  flight  conditions  to  support  a  test  by  Captain  John  Anselmo,  Summer  1990. 

C 

C 

COMMON /MAIN1 /CMACH , ALPHA , NFLWS , IFLAG, CNTR01D ( 2000, 3 ) , 

1  DIRECOS ( 2000 , 3) , AREA (2000) , CP (2000) 

C 

I FLAG* 1 
C 

C  *****  IDENTIFY  THE  SOURCE  OF  AERO  DATA  ***** 

C 

OPEN ( 2 , FILE* ' B1FUS2 . DMP ' , STATUS* 'OLD ' ) 

OPEN ( 3 , FILE* ' CANOPY . GEO ' , STATUS* ' NEW ' ) 

OPEN ( 4 , FILE* ' CANOPY . ARO ' , STATUS* ' NEW' ) 

1=0 

4010  READ  ( 2 , 4020 , END=4030 )  QU 
4020  FORMAT (A4) 

1*1+1 

IF  (QU.EQ. 'QUAD' )  GOTO  4050 
IF  (I.LT.2000)  GOTO  4010 
C 

C  *****  THE  AERO  SOURCE  CANNOT  BE  IDENTIFIED  ***** 

C 

4030  WRITE  (*,4040) 

4040  FORMAT ( '  FILE  2  was  not  identified  as  a  QUADPAN  dump  file') 

STOP  01 
C 


C  *****  THE  AERO  DATA  HAS  BEEN  IDENTIFIED  ***** 
C  *****  AS  A  QUADPAN  DUMP  FILE  ***** 

C 


4050  WRITE  (*,4060) 

4060  FORMAT ( '  FILE  2  has  been  identified  as  a  QUADPAN  file.') 
9000  CONTINUE 


C 

c  *****  read  CP  FOR  ONE  FLIGHT  CONDITION  ***** 
C 

CALL  READQUAD 
C 

C  *****  WRITE  CP  FOR  ONE  FLIGHT  CONDITION  ***** 
C 


CALL  WRITECP 
C 

C  *****  ip  LAST  FLOW  CONDITION,  THEN  STOP  PROGRAM  * 
C  *****  if  NOT,  THEN  CONTINUE  WITH  NEXT  FLOW  CONDITION 


G-l 


ooo  noonno  ooooo  oooo  oooo  oooo  oo 


IF  ( I FLAG  .EQ.  NFLWS )  STOP  02 

IFLAG-IFLAG+1 

GO  TO  9000 

END 


SUBROUTINE  READQUAD 

*****  THIS  SUBROUTINE  READS  A  QUADPAN  DUMP  FILE,  *’ 
*****  EXTRACTS  THE  REQUIRED  INFORMATION  ***** 

COMMON /MAIN1 /CMACH , ALPHA , NFLHS , I FLAG , CNTROID ( 2000 , 3 ) , 

1  DIRECOS (2000, 3 ) , AREA ( 2000 ) , CP (2000) 

*****  read  IN  GEOMETRY  DATA  ***** 

*****  ONLY  IF  THIS  IS  THE  FIRST  PASS  ***** 

IF  ( I FLAG  .GT.  1)  GO  TO  170 

*****  DETERMINE  NUMBER  OF  FLOW  CONDITIONS  AND  NUMBER  OF 
*****  ELEMENTS  IN  MODEL  ***** 

1  READ  (2,10)  GL 
10  FORMAT  (A4) 

IF  (GL  .NE.  'GLOB')  GO  TO  1 
READ  (2,20)  ITOTL,  NFLWS 
20  FORMAT  ( //2 ( 12X, 16) ) 


*****  DETERMINE  POSITIONS  OF  CONTROL  POINTS,  ***** 

*****  DIRECTION  COSINES,  AND  AREAS  OF  ELEMENTS  ***** 

70  READ  (2,10)  GE 

IF  (GE  .NE.  'GEOM' )  GO  TO  70 
READ  (2,90) 

90  FORMAT  (///) 

DO  110  K-l, ITOTL 

READ  (2,100)  ( CNTROID (K,J) , J-l, 3) ,  (DIRECOS (K,KK) ,KK-1, 3) ,  AREA(K) 
100  FORMAT  (/5G15.7,/ , 2G15.7) 

110  CONTINUE 

*****  THE  FOLLOWING  ARE  TO  BE  REPEATED  ***** 

*****  FOR  EACH  FLOW  CONDITION  ***** 

*****  DETERMINE  MACH  NUMBER  ***** 

170  READ  (2,10)  OU 

IF  (OU  .NE.  'OUTP' )  GO  TO  170 
READ  (2,180)  CMACH, ALPHA 
180  FORMAT  (/////2G15.7) 

*****  DETERMINE  COEFFICIENTS  OF  PRESSURE  ***** 

190  READ  (2,10)  PR 

IF  (PR  .NE,  'PRES')  GO  TO  190 
READ  (2,210) 

210  FORMAT  (//) 

DO  230  K-l, 800 
READ  (2,220)  CP(K) 

220  FORMAT  (//G15.7) 

230  CONTINUE 


G-2 


onoo 


RETURN 

END 


SUBROUTINE  WRITECP 

*****  THIS  SUBROUTINE  WILL  WRITE  GEOMETRY  Af’D  PRESSURE  ***** 
*****  COEFFICIENTS  ON  AN  OUTPUT  FILE  ***** 

DIMENSION  OUT ( 7 ) 

COMMON /MAI N 1 / CMACH , ALPHA , NFLWS  f I FLAG , CNTRO I D ( 2 000 , 3 ) , 

1  DIRECOS ( 2000 , 3 ) , AREA (2000 ) , CP ( 2000 ) 

IF  ( IFLAG.GT. 1 )  GO  TO  30 
WRITE (3, 1000) 

1000  FORMAT ('  GEOMETRY  B1  CANOPY'/ 

1  '  I , J, XBAR, YBAR, ZBAR,XCOS, YCOS, 2 COS, AREA' ) 

C 

DO  20  1*1,20 

DO  20  J*l, 40 

K*40*(I-1)+J 

DO  21  L*l,3 

OUT ( L ) =CNTROID ( K ,  L ) 

21  OUT ( L+3 ) =D IRECOS ( K , L ) 

OUT(7)=AREA(K) 

20  WRITE(3, 2000)1, J, OUT 
2000  FORMAT ( 214 , 3F10 . 4 , 3F10 . 7 , F10 . 6 ) 

C 

30  CONTINUE 

WRITE (4,3000) CMACH , ALPHA 

3000  FORMAT ( '  B1  CANOPY  CP  FOR  MACH  *'F6.3'  AND  ALPHA  *'F5.2) 

C 

DO  40  1*1,20 
WRITE (4, 4000) I 

4000  FORMAT ( '  CP(I,J)  FOR  I  -  '13'  AND  J  -  1  THROUGH  40') 
KFIRST-40* ( 1-1 ) +1 
KLAST*KFIRST+39 

40  WRITE(4, 5000) (CP(K) , K=KFIRST , KLAST ) 

5000  FORMAT (8F10. 6) 

RETURN 

END 


G-3 


